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CHAPTER  1 
THEORY 

1.1     Particle  Physics 

Particle  physics  is  the  branch  of  physics  dedicated  to  the  study  of  matter 
a.nd  energy  at  the  most  fundamental  level.  This  means  that  the  job  of  a  particle 
physicist  is  to  identify  the  smallest  constituents  of  matter  and  describe  how  they 
interact  with  each  other. 

Humankind  has  been  interested  in  this  subject  since  ancient  times.  Two  of  the 
first  particle  physicists  in  recorded  history  are  the  Greek  thinkers  Empedocles  and 
Democritus  from  the  fifth  century  BC.  Empedocles  stated  that  our  complex  world 
was  made  from  combining  four  fundamental  elements  (earth,  air,  fire  and  water)  in 
different  proportions.  Democritus  on  the  other  hand,  believed  that  the  apparently 
continuous  objects  in  the  natural  world  were  not  really  continuous,  but  made  from 
voids  and  indivisible  particles  called  atoms. 

Over  the  last  half  century  Particle  physics  has  advanced  tremendously  and  we 
now  have  a  beautiful,  but  incomplete,  theory  firmly  grounded  on  experiment  that 
describes  the  fundamental  constituents  of  matter  and  how  they  interact  with  each 
other.  This  theory  is  called  "The  Standard  Model." 
1.1.1     The  Standard  Model 

According  to  the  Standard  Model,  the  fundamental  building  blocks  of  matter 
are  point  like  particles  which  interact  with  each  other  in  as  many  as  three  different 
ways.  Each  type  of  interaction,  or  force,  is  itself  carried  by  point  like  particles 


called  force  carriers.  The  particles  which  are  force  carriers  are  bosons  1    and  are 
collectively  called  gauge  bosons  because  they  are  needed  for  the  theory  to  be  gauge 
invariant.  The  non-force  carrying  particles  are  fermions.  2 

The  three  interactions  described  by  the  Standard  Model  are  called  the 
electromagnetic  force,  the  weak  force,  and  the  strong  force.  The  force  of  gravity  is 
not  included  in  the  Standard  Model,  and  this  is  one  of  the  reasons  the  Standard 
Model  is  not  yet  complete.  The  electromagnetic  force  carriers  are  photons,  the 
weak  force  carriers  are  the  Z°,  W+,  and  W~  particles,  and  finally,  the  strong  force 
is  carried  by  eight  kinds  of  gluons.  Table  1-1  summarizes  the  situation. 
Table  1-1:  Gauge  bosons  and  the  force  they  carry. 


Symbol 

Name 

Force  Carried 

7 

Photon 

Electromagnetic 

Z° 

Z° 

Weak 

w+ 

w+ 

Weak 

w~ 

w~ 

Weak 

9 

gluon 

Strong 

Particles  which  interact  through  a  particular  force  are  said  to  couple  to  it 
and  to  carry  an  associated  charge.  The  nomenclature  is  as  follows,  particles  that 
interact  through  the  electromagnetic  force  have  an  electromagnetic  charge  called 
electric  charge,  those  which  interact  through  the  weak  force  have  a  weak  charge 
called  weak  isospin,  and  those  which  interact  strongly  have  a  strong  charge  called 


color. 


The  rest  of  the  particles  in  the  Standard  Model  which  are  not  force  carriers, 
the  fermions,  are  subdivided  into  different  groups  depending  on  their  properties 
(See  table  1-2)  ,  reminiscent  of  the  way  chemists  organized  the  elements  into 


1  Bosons  are  defined  as  particles  with  integer  spin  in  quantum  mechanics. 

2  Fermions  are  defined  as  particles  with  half  odd  fractional  spin. 


the  periodic  table  during  the  second  half  of  the  nineteenth  century.  Fermions  are 

divided  into  quarks  (generically  represented  by  the  symbol  q)  and  leptons.  The 

main  difference  between  these  two  groups  is  that  quarks  interact  through  the  strong 

force  while  leptons  do  not.  There  are  six  types  of  leptons  and  six  types  of  quarks 

(also  called  the  six  quark  flavors)  which  are  grouped  into  three  generations.  Each 

generation  consist  of  two  quarks  and  two  leptons.  All  three  generations  replicate 

the  same  set  of  force  charges,  the  main  difference  between  generations  is  the  mass 

of  the  particles  (for  example,  the  ratio  of  the  masses  of  e  :  /i  :  T  is  1  :  200  :  3500). 

For  each  fermion  there  is  an  anti-fermion  with  equal  mass  and  spin  and  opposite 

charge. 

Table  1-2:  Fermion  symbols  classified  into  quarks,  leptons,  and  the  three  genera- 
tions along  with  the  generational  common  charges. 


First 
Generation 

Second 
Generation 

Third 
Generation 

Electric 
Charge 

Weak 
Charge 

Has 
Color 

Quarks 

u 
d 

c 

s 

t 

b 

+2/3 
-1/3 

+1/2 
-1/2 

Yes 
Yes 

Leptons 

e 

T 

0 

-1 

+1/2 
-1/2 

No 

No 

1.1.2     Quantum  Chromodynamics 

Quantum  chromodynamics  (QCD)  is  the  part  of  the  Standard  Model  that 
describes  the  strong  force.  QCD  is  based  on  local  gauge  invariance  and  color 
symmetry.  There  are  three  possible  color  charges  for  quarks  called  r  (red),  b  (blue), 
and  g  (green).  Anti-quarks  have  opposite  colors  called  f,b,g.  The  strong  force 
between  quarks  only  depends  on  their  colors  and  is  independent  of  their  flavor. 

A  very  important  characteristic  of  the  strong  force  is  that  gluons  themselves 
carry  a  color  charge  and  an  anti-color  charge  so  they  can  interact  with  other  quarks 
through  the  strong  force  and  change  their  color.  Since  there  are  three  colors  and 
three  anti-colors  one  might  think  that  there  are  nine  different  gluons.  However, 
there  is  one  linear  combination  of  color  anti-color  states  that  has  no  net  color  and 


leaves  a  quark  unchanged.  There  are  therefore  9-1  =  8  gluons.  The  8  individual 
gluon  color  states  can  be  written  as  follows, 

|1>  =  -={rb  +  bf) 
v2 

12  >  =  ^{rb-bf) 

|3>  =  -=[rf-bb) 
V2 

|4>=  -^{rg  +  gf) 

V2  (i.i) 

|5>=  --Urg-gf) 

|6>=  -f={bg  +  gb) 

V2 

|7  >  =  ~(bg  -  gb) 

\S>  =  -={rf  +  gg-2bb) 

vo 

and  represent  the  eight  different  gluons  that  exist  in  nature.  The  single  and  unique 
color  state  left  out  is  called  a  color  singlet, 

|9>=  -^(^  +  66  +  gg)  (1-2) 

v3 

which  is  invariant  under  a  redefinition  of  the  color  (a  rotation  in  color  space).  In 
group  theory  this  decomposition  of  the  color  states  into  an  octet  and  a  singlet  is 
denoted  by  3  <g>  3  =  8  ©  1.  It  is  worth  noting  here  that  a  colorless  sate,  such  as  |3  > 
or  1 8  >,  is  not  necessarily  a  color  singlet  state. 

This  situation  is  analogous  to  the  perhaps  more  familiar  example  of  two  spin 
1/2  particles.  Each  particle  can  have  their  spin  up  (t)  or  down  (|)  corresponding  to 
four  possible  combinations,  each  giving  a  total  spin  S  =  0  or  1.  The  S  =  1  states 


form  a  triplet, 


|1>  =  |TT> 


i 


|2  >  =  -^(|  U>  +|  it>) 


(1.3) 


|3  >  =  |  U> 
and  there  is  a  singlet  state  with  S  =  0, 

1 


|4>=^=(|U>H1T>). 


(1.4) 


Since  gluons  themselves  carry  color  they  can  interact  with  each  other  through 
the  strong  force.  This  interaction  among  the  force  carriers  is  unique  to  the  strong 
force  and,  when  included  in  perturbative  QCD  calculations,  leads  to  two  important 
properties  observed  in  nature  called  "asymptotic  freedom"  and  "color  confinement." 
Asymptotic  freedom  means  that  the  interaction  gets  weaker  at  short  distances. 
Color  confinement  is  the  requirement  that  observed  states  have  neutral  color,  or  in 
other  words,  they  must  be  in  a  color  singlet  state. 

Confinement  explains  why  free  quarks  or  free  gluons,  which  have  a  net  color 
charge,  have  never  been  observed.  It  also  explains  why  no  fractional  charged 
particles  made  from  a  qq  bound  state  have  never  been  observed  since  it  is  not 
possible  to  construct  a  color  singlet  for  such  a  state  (in  group  theory  terms 
3  <g>  3  =  6  ©  3  where  we  have  a  sextet  and  a  triplet,  but  no  singlet).  On  the  other 
hand,  color  singlets  can  be  constructed  for  a  qq  or  qqq  system.  The  color  singlet  for 
qq  is  simply  the  state  shown  on  the  right  side  of  Equation  1.2.  The  color  singlet  for 
qqq  can  be  obtained  from  the  decomposition  3®3®3--=10©8©8©l  and  is, 


\qqq  > color  sin9iet=  -j=  {rgb  -  grb  +  brg  -  bgr  +  gbr  -  rbg).  (1.5) 

Particles  that  are  bound  states  of  qq  are  abundant  in  nature  and  are  called 
mesons,  those  that  are  bound  states  of  qqq  are  called  baryons  and  also  abound  in 


nature.  Both  groups  axe  collectively  called  hadrons.  There  have  been  hundreds 
of  different  hadrons  observed,  confirming  the  validity  of  QCD  and  the  quark 
flavors.  Nevertheless,  QCD  leaves  room  for  more  possibilities.  A  bound  state  of  two 
gluons  gg  (sometimes  called  a  glueball)  can  be  in  a  color  singlet  and  in  principle 
could  be  observed.  Although  glueballs  are  allowed  by  QCD  there  is  no  convincing 
experimental  observation  of  one.  Another  possibility  are  bound  states  that  are 
a  mixture  of  the  previous  states,  such  as  aqq  +  (3gg  with  arbitrary  a,  (5.  These 
states  are  called  "hybrid  mesons"  or  simply  "hybrids".  It  is  believed  [1],  [2],  [3]  that 
hybrids  are  necessary  to  explain  the  spectrum  of  light  mesons  between  1-2   GeV. 

As  with  any  quantum  mechanical  theory,  the  spectrum  of  bound  states  is 
a  fundamental  test. 3    Glueballs  are  allowed  by  QCD,  yet  there  is  no  conclu- 
sive experimental  evidence  of  their  observation,  despite  intense  experimental 
searches  [4],  [5],  [6]  complemented  by  lattice  QCD  calculations  [7],  [8]  and  other 
theoretical  contributions  like  bag  models  [9],  flux-tube  models  [10],  QCD  sum 
rules  [11],  weakly  bound  bound-state  models  [12],  and  QCD  factorization  formalism 
models  [13].  Physicists  cannot  be  sure  they  understand  QCD  until  such  states  are 
observed,  or  until  they  can  explain  why  we  cannot  observe  them. 
1.1.3     Introduction  to  the  Radiative  Decays  of  the  T(1S) 

The  T(1S)  is  a  meson  composed  of  bb  quarks.  Mesons  of  this  kind,  composed 
of  a  quark  and  an  anti-quark  of  the  same  flavor,  are  in  general  called  quarkonia. 

There  is  a  convention  behind  the  name  of  the  T(1S)  .  The  T  symbol  is 
reserved  for  particles  composed  of  bb  where  the  combined  spin  of  the  quark  and 
anti-quark  is  1.  The  "IS"  symbol  is  borrowed  from  atomic  spectroscopy  with  the 
"1"  meaning  that  the  bb  pair  are  in  the  lowest-energy  bound-state,  and  the  "S" 


3  A  good  example  of  this  is  the  successful  description  of  the  hydrogen  atom's 
spectrum  by  quantum  mechanics. 


meaning  that  the  bb  have  a  relative  angular  momentum  L  ==  0.  A  description  of 
particle  naming  conventions  can  be  found  in  [14]. 

The  T(1S)  is  unstable,  existing  for  only  about  10-24  seconds  after  which  it 
decays  into  daughter  particles  (which  in  turn  decay  themselves  if  they  are  not 
stable).  The  term  "radiative  decay"  is  reserved  to  any  T(1S)  decay  where  one  of 
the  stable  daughters  is  a  photon. 

The  different  ways  the  T(1S)  can  decay  must  obey  the  symmetries  in  nature. 
Examples  of  such  symmetries  are  "parity"  and  "charge  conjugation" ,  both  of  which 
will  be  described  soon.  Symmetries  are  very  important  in  the  standard  model.  The 
usefulness  of  symmetries  can  be  seen  in  Noethers'  theorem,  which  states  that  for 
every  symmetry  there  is  an  associated  conserved  quantity.  Examples  of  Noethers' 
theorem  are  the  conservation  of  momentum  connected  to  translational  invariance 
and  the  conservation  of  angular  momentum  associated  with  rotational  invariance. 

In  the  next  sections  we  will  use  the  parity  and  charge  conjugation  symmetries 
to  find  T(1S)  decays  which  conserve  the  associated  symmetry  constants  and  are 
allowed  by  nature.  In  particular,  we  will  show  that  the  radiative  decay  of  the 
T(1S)  through  a  photon  and  two  gluons  (see  Figure  1-1)  is  allowed.  The  key 
observation  is  that  the  two  gluons  must  be  in  a  color  singlet  since  both  the  T(1S) 
and  the  radiated  photon  have  no  color  and  color  must  be  conserved.  This  means 
that  the  two  gluons  satisfy  color  confinement  and  could  form  a  glueball,  although 
more  conventional  meson  states,  or  hadrons  in  no  bound  state  at  all,  are  also 
possible  outcomes.  Regardless  of  what  the  gluons  do,  their  energy  will  eventually 
manifests  itself  as  hadrons.  Sometimes,  two  charged  hadrons  of  opposite  charge  will 
emerge.  This  work  examines  those  two  hadrons  from  a  radiative  T(1S)  decay  to 
experimentally  probe  the  gg  spectrum. 


Parity 

The  parity  operator,  P,  reverses  the  sign  of  an  object's  spatial  coordinates. 
Consider  a  particle  \a  >  with  a  wave  function  ^a(x,t).  By  the  definition  of  the 

parity  operator, 

P<l>(x,t)  =  Pa<S>a(-x,t)  (1-6) 

where  Pa  is  a  constant  phase  factor.  If  we  consider  an  eigenfunction  of  momentum 


M^t) 


,i(p-x-Et) 


(1.7) 


then 

F9p{x,  t)  =  Pa^p(-x,  t)  -  P.f  Jfc  t),  (1.8) 

so  that  any  particle  at  rest,  with  p*0,  remains  unchanged  up  to  a  multiplicative 
number,  Pa,  under  the  parity  operator.  States  with  this  property  are  called 
eigenstates  with  eigenvalue  Pa.  Pa  is  also  called  the  intrinsic  parity  of  particle  a, 
or  more  usually  just  the  parity  of  particle  a,  with  the  words  at  rest  left  implicit. 
Since  two  successive  parity  transformations  leave  the  system  unchanged,  P*  «=  1, 
implying  that  the  possible  values  for  the  parity  eigenvalue  are  Pa  =  ±1. 

In  addition  to  a  particle  at  rest,  a  particle  with  definite  orbital  angular 
momentum  is  also  en  eigenstate  of  parity.  The  wave  function  for  such  a  particle  in 
spherical  coordinates  is, 


^nlm(x,t)  =  Rnl(r)Yr(e,<P), 


[1.9) 


where  (r,  6,  (j))  are  spherical  polar  coordinates,  Rni{r)  is  a  function  of  the  radial 
variable  r  only,  and  the  Ylm(0,4>)  is  a  spherical  harmonic. 

The  spherical  harmonics  are  well  known  functions  which  have  the  following 
property, 

Yne,  4)  =  (-i)f*r(*  -  m + <i>)-  (L1°) 


Hence 

PVnlm(Z,t)  =  Pa^>nlm(-X,t)  =  P^-l^W*  ,*)  (!•") 

proving  that  a  particle  with  a  definite  orbital  angular  momentum  f  is  indeed  an 
eigenstate  of  the  parity  operator  with  eigenvalue  Pa(-if. 

The  parities  of  the  fundamental  fermions  cannot  be  measured  or  derived.  All 
that  nature  requires  is  that  the  parity  of  a  fermion  be  opposite  to  that  of  an  anti- 
fermion.  As  a  matter  of  convention  fermions  are  assigned  P  =  +1  and  anti-fermions 
are  assigned  P  =  -1.  In  contrast,  the  parities  of  the  photon  and  gluon  can  be 
derived  by  applying  P  to  the  field  equations  resulting  in  P7  =  — 1  and  Pg  =  —  1. 

The  T(1S)  has  P  =  PbP-b(-l)L  =  -1  since  L  =  0. 

Parity  is  a  good  quantum  number  because  it  is  a  symmetry  of  the  strong  and 
electromagnetic  force.  This  means,  that  in  any  reaction  involving  these  forces, 
parity  must  be  conserved. 
Charge  Conjugation 

Charge  conjugation  is  simply  the  operation  which  replaces  all  particles  by 
their  anti-particles.  In  quantum  mechanics  the  charge  conjugation  operator  is 
represented  by  C.  For  any  particle  \a  >  we  can  write, 


C\a  >=  ca\a  > 


[1.12) 


where  ca  is  a  phase  factor.  If  we  let  the  C  operator  act  twice  to  recover  the  original 
state  \a  >, 

\a  >=  C2\a  >=  C(ca\a  >)  =  caC\a  >=  cac-a\a  >  (1.13) 

which  shows  that  caCa  =  1.  If  (and  only  if)  a  is  its  own  anti-particle,  it  is  an 
eigenstate  of  C.  The  possible  eigenvalues  are  limited  to  C  —  ca  —  c&  =  ±1. 

All  systems  composed  of  a  the  same  fermion  and  an  anti-fermion  pair  are 
eigenstates  of  C  with  eigenvalue  C  =  (-l)(L+5).  This  factor  can  be  understood 
because  of  the  need  to  exchange  both  particles'  position  and  spin  to  recover  the 
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original  state  after  the  charge  conjugation  operator  is  applied.  Exchanging  the 
particles'  position  gives  a  factor  of  (-1)L  as  was  shown  in  the  previous  section, 
exchanging  the  particles  spin  gives  a  factor  of  (-l)s+1  as  can  be  verified  by 
inspecting  Equations  1.3  and  1.4,  and  a  factor  of  (-1)  which  arises  in  quantum  field 
theory  whenever  fermions  and  anti-fermions  are  interchanged.  With  this  result  we 
can  calculate  the  charge  conjugation  eigenvalue  for  the  T(1S)  and  obtain  C  —  -1 
since  L  +  S  —  1. 

The  photon  is  an  eigenstate  of  C  since  it  is  its  own  anti-particle.  The  C 
eigenvalue  for  the  photon  can  be  derived  by  inserting  C  into  the  field  equations  and 
is  C7  =  -1. 

Finally,  we  consider  a  system  composed  of  two  gluons  that  are  in  a  color 
singlet.  The  two  gluons  are  bosons  and  they  must  have  a  symmetric  wave  function, 
^G  under  a  gx  *-*  g2  exchange.  Under  this  exchange,  the  orbital  angular  momentum 
part  of  the  wave  function  contributes  with  a  factor  (-1)L,  the  spin  part  of  the 
wave  function  for  two  spin  1  particles  contributes  with  a  factor  of  (-l)s,  and  the 
color  singlet  part  of  the  wave  function  contributes  with  a  factor  of  +1  since  it  is 
symmetric.  To  ensure  that  $G  is  symmetric  we  need  L  +  S  to  be  even.  This  implies 
that  C  =  (-1)L+S  =  +1.  The  L  =  0  and  L  =  1  possible  gg  bound  states  are  shown 
in  Table  1-3.  The  JPC  =  1~+  is  peculiar  because  it  is  impossible  for  a  qq  system  to 
have  these  quantum  numbers.  If  this  state  is  ever  observed,  it  must  be  a  glueball. 
Experimental  searches  for  such  a  state  have  been  done,  for  example,  in  [15]. 

Table  1-3:  Possible  gg  bound  states  with  L  «  0  or  L  =  1.  The  possible  quantum 
numbers  are  limited  by  the  condition  that  J  =  L  +  S  be  even,  which  is  needed  to 
ensure  a  symmetric  wave  function  for  the  two  gauge  bosons. 


L 

S 

jPC 

0 
0 
1 

0 
2 
1 

0++ 

2++ 
0-+,l-+,2-+ 
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Charge  conjugation  is  a  symmetry  of  the  strong  and  electromagnetic  force.  For 
those  particles  that  are  eigenstates  of  C,  C  is  a  good  quantum  number  because  in 
any  reaction  involving  these  forces  C  must  be  conserved. 

The  P  and  C  values  of  various  particles  used  in  this  analysis  are  shown  in 
Table  1-4  along  with  their  quark  composition,  orbital  angular  momentum,  and 

internal  spin. 

Table  1-4:  Symbol,  name,  quark  composition,  angular  momentum  (L),  internal  spin 
(S),  parity  (P),  and  charge  conjugation  eigenvalues  (C)  for  a  few  of  the  particles 
used  in  this  analysis. 


Symbol 

Name 

Quark  Composition 

L 

S 

P 

c 

T(1S) 

Upsilon(lS) 

66 

0 

1 

-1 

-1 

7!-+ 

Pion 

ud 

0 

0 

-1 

X 

7T~ 

Pion 

du 

0 

0 

-1 

X 

K+ 

Kaon 

us 

0 

0 

-1 

X 

K- 

Kaon 

su 

0 

0 

-1 

X 

P 

Proton 

uud 

0 

1/2 

+1 

X 

V 

Anti-proton 

i                 uud 

0 

1/2 

-1 

X 

1 

Photon 

X 

X 

1 

-1 

-1 

Possible  Decays  of  the  T(1S) 

At  this  point  we  can  understand  the  different  possible  ways  the  T(1S)  can 
decay.  The  Possible  decays  are  limited  because  the  strong  and  electromagnetic 
force  must  conserve  color,  C  and  P. 

The  simplest  possibility  is  for  the  66  pair  to  interact  electro-magnetically  and 
annihilate  into  one  virtual  photon  4  .  This  is  allowed  by  parity,  charge  conjugation 
symmetry  and  color  conservation.  The  decay  of  the  T(1S)  to  one  gluon  is  not 
allowed  by  color  conservation  and  is  therefore  forbidden.  T(1S)  decays  to  two 
photons  are  forbidden  by  charge  conjugation.  T(1S)  decays  to  two  gluons  are 


4  Such  a  photon  is  called  virtual  because  it  cannot  conserve  the  4-momentum  of 
66  and  is  unstable,  only  existing  for  a  brief  period  of  time,  as  allowed  by  the  uncer- 
tainty principle,  after  which  it  decays. 
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also  forbidden  by  charge  conjugation.  T(1S)  decays  to  3  gluons  are  allowed  (3 
gluons  can  form  a  color  singlet).  T(1S)  decays  to  3  photons  are  also  allowed,  but 
are  largely  suppressed  by  decays  to  one  photon  since  3  successive  electromagnetic 
interactions  are  much  less  likely  to  occur  than  a  single  one.  Finally,  T(1S)  decays 
two  one  photon  and  two  gluons  are  also  allowed  under  the  condition  that  the  two 
gluons  be  in  a  color  singlet  state. 

The  three  different  possible  T(1S)  decays  with  least  amount  of  interactions 
(also  called  lowest  order  decays)  are  shown  in  Figure  1-1. 


(a) 


(b) 


(c) 


Figure  1-1:  Lowest  order  decays  of  the  T(1S)  allowed  by  color  conservation,  charge 
conjugation  symmetry,  and  parity,  (a)  Shows  the  decay  into  three  gluons,  (b)  shows 
a  radiative  decay,  and  (c)  shows  the  electromagnetic  decay  through  a  virtual  pho- 
ton that  in  turn  decays  electromagnetically  into  a  pair  of  charged  fundamental 
particles,  such  as  quarks  or  charged  leptons  (the  charged  leptons  are  represented  by 
the  symbol  I). 


Observable  Resonances 

In  this  work  we  search  for  a  resonance  X  produced  in  a  radiative  T(1S)  decay, 


T(1S)  -+  iX. 


(-1-14) 
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Using  charge  conjugation  symmetry  on  both  sides, 

-lL-lCx.  (1.15) 

Therefore,  Cx  -  +1. 

In  order  for  us  to  observe  X  it  must  decay  into  two  charged  hadrons, 

X  -»  h+h~,  (1.16) 


where  ft,  =  tt,  JC,  p  are  the  hadrons  whose  momentum  we  are  going  to  measure. 
Applying  charge  conjugation  to  this  last  decay, 


+1  =  (-1) 


L+S 


(1.17) 


where  L  and  S  are  respectively  the  angular  momentum  and  spin  of  the  h+  h~ 
system.  This  last  equation  implies  that  L  +  S  must  be  even.  This  has  consequences 
for  the  possible  X  parities  we  can  observe.  By  parity  conservation  in  1.16, 


P 


x  -- 


■1)L(-1)(-1)  -  (-1)J 


(1.18) 


For  h  =  7r,  K:  S  =  0  and  L  ~  J  must  be  even,  which  implies  that  Px  =  -Hi.  For 
h  =  p,  S  can  be  0  or  1  and  particles  with  both  positive  and  negative  parities  can  be 

detected. 

Table  1-5:  Possible  5,  L,  J,  P  and  C  values  for  X  from  the  radiative  decay 
T(IS)  — >  7X  reconstructed  in  different  decay  modes. 


Decay  Mode 

s 

L 

J 

P 

C 

X  — >  7T+7T~ 

0 

even 

even 

+1 

+1 

x  ->  /c+x- 

0 

even 

even 

+1 

+1 

X  —+  pp 

0 

even 

even 

■fl 

+1 

X^'pp 

1 

odd 

even  and  odd 

-1 

+1 
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1.2     Radiative  Decays  of  Quarkonia  Overview 

Theoretical  models  exist  for  glueball  production  in  quarkonia  decay  [12] 
and  for  the  glueball  spectrum.  For  example,  a  quenched  lattice  calculation  [7] 
predicts  a  Jpc  =  24+  glueball 5   in  the  2.2  GeV/c2  mass  region  (see  Figure  1-2). 
According  to  Table  1-5,  in  the  charged  pion  and  kaon  modes,  we  are  limited  to 
detect  glueballs  in  the  leftmost  column  where  P  =  C  ==  +1  of  Figure  1-2,  while  in 
the  proton  mode  we  are  restricted  to  the  two  left  most  columns. 
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Figure  1-2:  Quenched  lattice  calculation  result  for  the  glueball  spectrum  for  differ- 
ent P  and  C  values.  The  mass  scale  is  shown  in  terms  of  a  scaling  parameter  from 
the  QCD  lattice  calculation  named  r0  on  the  left  and  in  absolute  terms  on  the  right 
by  taking  r^1  =  410  MeV. 


5  Here  J  stands  for  the  internal  angular  momentum  (spin)  of  the  glueball, 

J  =  L  +  S. 
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For  some  time  bound  states  not  considered  to  be  pure  glueballs,  such  as 
/2(1270),  /4(2050),  r),  and  rf  have  been  observed  in  J/ip  radiative  decays  at 
the  10"3  production  level  [16]  6  .  In  1996,  the  BBS  collaboration  claimed  the 
observation  of  a  resonance  they  called  the  /,/(2220)  particle  in  the  radiative  J  ftp 
system  at  the  1(T5  level  [17]  (see  Figure  1-3).  A  lot  of  excitement  was  generated  at 
the  time  because  it  is  possible  to  interpret  the  /j(2220)  as  a  glueball.  However,  this 
result  has  not  been  confirmed. 

27/0301-004 

i 1 1 1 1 r~- 

(a)  x*ri~    - 


^VcP-c 


2.2  2.4 

Mass  (GeV) 


Figure  1-3:  The  mass  spectrum  obtained  by  the  BES  collaboration  in  radiative 
J/ijj  decays  into  different  hadronic  modes. 


Key  to  identifying  a  particle  as  a  glueball  are  (a)  suppressed  production  in 
two-photon  collisions  (unlike  quarks,  gluons  don't  carry  electric  charge  and  do  not 


6  The  J/ip  particle  is  a  cc  bound  state.  Since  the  strong  force  is  flavor  blind  the 
situation  in  J/ip  decays  is  in  principle  similar  to  that  in  T(1S)  decays. 
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couple  to  photons),  and  (b)  flavor  symmetric  decays,  since  a  pure  glueball  has  no 
valence  quarks. 

CLEO  has  already  done  several  studies  of  radiative  decays  of  the  T(lS'). 
Naively,  one  expects  these  types  of  decays  to  be  suppressed  by  a  factor 


\(qb/qc)(rrijmb)}2  m  0.025 


(1.19) 


with  respect  to  J/tp  radiative  decays.  Phis  comes  from  noticing  that  the  quark- 
photon  coupling  is  proportional  to  the  electric  charge  and  the  quark  propagator 
is  approximately  1/m  for  low  momentum  quarks.  In  1999,  CLEO  made  the  first 
observation  of  a  two-body  T(IS')  radiative  decay  [18].  The  spin  of  the  observed 
resonance  could  not  be  measured,  but  its  mass  and  width  where  consistent  with  the 
/2(1270)  particle.  Under  this  assumption,  comparing  the  measured  branching  ratio 
of  T(1S)  -►  7/2(2220)  to  the  measured  branching  fraction  of  the  J/ip  -*  7/2  (1270), 
a  suppression  factor  of  0.06  ±  0.03  was  obtained.  After  the  BES  result  for  the 
f 'j  (2220)  in  radiative  J/ip  decays,  a  corresponding  search  was  performed  by  CLEO 
in  the  radiative  T(15)  system  [19].  This  analysis  put  limits  on  the  /j(2220) 
production  in  radiative  T(1S)  decays. 

In  this  work  we  are  privileged  to  have  available  the  largest  collection  of 
radiative  T(1S)  decays  in  the  world.  With  it,  we  can  study  the  structure  of  color 
singlet  gg  hadronization,  and  shed  more  light  on  the  /j(2220)  result  from  BES. 


CHAPTER  2 
EXPERIMENTAL  APPARATUS 

To  carry  out  our  study  of  the  di-hadron  spectrum  we  need  to  first  produce 
the  T(1S)  resonance  and  secondly  observe  its  daughter  particles  flying  away  at 
relativists  speeds.  These  two  tasks  axe  respectively  accomplished  by  the  Cornell 
Electron  Storage  Ring  and  the  CLEO  III  detector. 

2.1     The  Cornell  Electron  Storage  Ring 
The  Cornell  Electron  Storage  Ring  (CESR),  located  at  Cornell  University,  is  a 
circular  particle  accelerator  that  produces  e+e~  collisions. 

In  order  to  produce  such  collisions  electrons  and  positrons  need  to  be  created, 
accelerated  and  stored.  CESR's  different  components,  shown  in  Figure  2-1,  have 
been  carrying  out  this  task  since  1979. 

A  typical  CESR  rim  begins  at  the  linear  accelerator  (LINAC)  where  electrons 
and  positrons  are  produced.  To  create  positrons,  electrons  are  evaporated  off  a 
filament  and  linearly  accelerated  by  electromagnetic  fields  towards  a  tungsten 
target.  The  collision  creates  a  spray  of  electrons,  positrons  and  photons.  The 
electrons  are  cleared  away  with  magnetic  fields  and  the  positrons  are  introduced 
into  the  synchrotron.  The  filling  procedure  is  identical,  except  that  the  tungsten 
target  is  removed. 

Once  the  electron  and  positron  beams  are  introduced  into  the  synchrotron, 
they  are  accelerated  to  the  operating  energy.  In  hour  case  they  are  accelerated  to 
the  point  where  their  combined  energy  is  the  T(1S)  mass,  9.46   GeV. 

Once  the  beams  are  at  the  desired  energy  they  are  transfered  to  the  storage 
ring,  where  they  will  remain  for  about  an  our.  At  one  point  in  the  storage  ring  the 
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beams  are  forced  to  croas  paths.  This  is  the  point  where  e+e    collisions  occur  1 
and  where  the  center  of  the  CLEO  III  detector  is  located.  If  the  accelerator  is 


•   BUNCH  OF  POSITRONS 
O  BUNCH  OF  ELECTRONS 


Figure  2-1:  The  Wilson  Laboratory  accelerator  located  about  40  feet  beneath 
Cornell  University's  Alumni  Fields. 

performing  well  a  high  collision  rate  results.  A  high  collision  rate  is  crucial  for  the 
success  of  an  accelerator  and  the  experiments  it  serves.  The  important  figure  is  the 
number  of  possible  collisions  per  second  per  unit  area;  this  is  called  the  luminosity. 
In  order  to  maximize  the  luminosity,  the  beams  are  focused  as  small  as  possible  at 


1  This  point  is  not  fixed  in  space,  but  varies  from  event  to  event  inside  a  small 
volume  of  space  called  the  interaction  region  (IR). 
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the  IR.  During  the  CLEO  III  installation  several  magnetic  quadruples  were  added 
to  CESR  to  improve  the  beam  focus.  CESR  has  consistently  outdone  itself  while 
collecting  luminosity  over  the  years  (see  Figure  2-2). 


2231001-009 

I r — 


II III  II 1 1 II I  ill  I  i 

Figure  2-2:  CESR  yearly  luminosity.  The  gaps  in  1995  and  1999  correspond  to 
down  times  when  the  CLEO  II.V  and  CLEO  III  detectors  where  being  installed. 

The  CLEO  III  detector  measures  the  time  integrated  luminosity  over  a  period 
of  time  by  counting  how  many  times  a  benchmark  process  occurs  during  a  certain 
time  interval  at  the  IR.  For  redundancy,  there  are  two  benchmark  processes  that 
are  used,  one  where  the  e+e~  particles  interact  to  produce  a  new  e+e -  pair, 
and  the  other  one  where  the  e+e~  annihilate  and  produce  two  photons.  Using 
the  known  cross-section  for  each  process,  the  number  of  events  is  converted  to  a 
luminosity  called  the  Bhabha  integrated  luminosity  for  the  first  process,  and  the  77 
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integrated  luminosity  for  the  second  one.  The  terra  "integrated"  is  sometimes  left 
out  and  the  total  luminosity  is  referred  to  as  simply  the  Bhabha  or  77  luminosity, 
with  the  time  integration  left  implicit. 

2.2     The  CLEO  III  Detector 

When  the  e+e~  collision  occurs,  the  two  particles  are  annihilated  we  enter 
the  world  of  particle  physics.  Nature  decides  what  to  do  with  the  energy  from 
the  annihilation.  We  have  no  chance  of  directly  observing  what  is  happening  at 
the  annihilation  point,  but  eventually  long  lived  semi-stable  particles  are  created 
that  fly  off  at  relativistic  speeds.  These  particles  carry  information  about  what 
happened  after  the  e+e~  collision,  and  can  tell  us  what  nature  did.  The  CLEO  III 
detector  has  the  important  mission  of  detecting  and  measuring  such  particles. 

As  one  can  see  in  Figure  2-3,  the  CLEO  III  detector  is  a  composite  of  many 
detector  elements.  These  sub-detectors  are  typically  arranged  as  concentric  cylin- 
ders. The  entire  detector  is  approximately  cube  shaped,  with  one  side  measuring 
about  6  meters,  and  weighs  over  500  thousand  kilograms. 
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Figure  2-3:  The  CLEO  III  detector. 
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As  implied  by  the  name,  CLEO  III  is  not  the  only  CLEO  detector.  CLEO 
III  was  preceded  by  CLEO  II.V,  CLEO^JI,  CLEO  I.V  and  the  original  CLEO 
detectors.  The  CLEO  III  detector  was  a  major  upgrade  compared  with  the 
previous  version  of  CLEO  [20],  [21],  and.  has  an  improved  particle  identification 
system  together  with  a  new  drift,  chamber  and  a  new  silicon  vertex  detector. 
2.2.1     Superconducting  Coil 

All  the  detector  subsystems  except  for  the  muon  chambers  are  located  inside 
a  superconduction  coil.  The  coil  remains  unchanged  since  CLEO  II.  It  is  kept  in 
a  superconducting  state  by  liquid  helium.  The  purpose  of  the  coil  is  to  provide  a 
1.5  Tesla  magnetic  field,  which  is  uniform  to  0.2%,  to  bend  the  paths  of  charged 
particles  in  the  detector.  By  measuring  how  much  much  a  charged  particle  bends, 
experimenters  can  measure  the  momentum  of  the  particle. 

The  coil  inner  radius  is  1.45m  and  its  outer  radius  is  1.55  m,  with  a  radial 
thickness  of  0.10  m.  The  total  length  of  the  coil  in  z  is  3.50  m.  It  is  wound  from 
a  5mm  x  16  mm  superconducting  cable  (Al  surrounding  Cu-NbTi  strands).  It 
is  wound  in  2  layers,  with  650  turns  per  layer,  on  an  aluminum  shell.  When  in 
operation  a  current  of  3300  amps  flows  through  the  coil. 
2.2.2     Tracking  System 

After  particles  from  the  interaction  point  pass  through  the  beam  pipe,  they 
begin  to  encounter  the  active  detector  elements  of  the  tracking  system.  There  are 
two  sub-detectors  responsible  for  tracking  the  curving  path  of  charged  particles. 
The  first  one  encountered  by  particles  is  the  silicon  vertex  detector,  and  the  second 
one  is  the  central  drift  chamber.  The  CLEO  III  tracking  system  is  responsible 
for  tracking  a  charged  particle's  path  and  measuring  its  momentum.  Typical 
momentum  resolution  is  0.3%  for  1  GeV  tracks.  The  tracking  system  also  measures 
ionization  energy  losses  with  an  accuracy  of  about  6%. 
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Silicon  Vertex  Detector 

The  silicon  vertex  detector  in  CLEO  III  [22],  also  called  SVD  III,  is  a  four- 
layer  barrel-only  structure  with  no  endcaps  that  surrounds  the  beam  pipe.  This 
detector  (see  Figure  2-4)  provides  four  (j>  and  four  z  measurements  covering  over 
over  93%  of  the  solid  angle.  The  average  radius  of  inner  surface  of  the  four  layers  is 
25  mm,  37.5  mm,  72  mm,  102  mm.  The  detector  is  constructed  from  447  identical 
double-sided  silicon  wafers,  each  27.0  nun  in  <£,  52.6  mm  in  z  and  0.3  nun  thick. 
The  wafers  are  instrumented  and  read  out  on  both  sides.  The  instrumentation  on 
each  side  consists  of  an  array  of  aluminum  strips  on  the  wafer  surface.  These  strips 
are  connected  to  preamplifiers  at  the  end  of  the  detector.  The  inner  side  has  512 
strips  in  the  z  direction  and  the  outer  side  has  512  in  the  cj>  direction.  Therefore 
each  wafer  contains  512+512  sensors.  1'he  447  wafers  are  arranged  in  the  4  layers, 
as  follows:  7  sections  in  0,  each  with  3  wafers  in  z,  total  ^  21  wafers  in  the  first 
layer;  10  sections  in  <£,  each  with  4  wafers  in  z,  total  =  40  wafers  in  the  second 
layer;  18  sections  in  4>.  each  with  7  wafers  in  z,  total  =  126  wafers  in  the  third 
layer;  26  sections  in  phi,  each  with  10  in  Z,  total  =  260  wafers  in  the  fourth  layer. 

Charged  particles  traversing  the  wafer  lose  energy  and  create  electron  hole 
pairs.  Approximately  3.6   eVis  required  to  create  a  single  electron-hole  pair.  The 
electrons  and  holes  then  travel  in  opposite  directions  in  the  electric  field  applied 
to  the  surfaces  of  the  wafers  until  they  end  up  on  the  aluminum  strips,  and  the 
detector  registers  a  "hit" .  When  combined  together,  the  hit  on  the  inner  side  of  a 
wafer  and  the  hit  on  the  outer  side  give  a  measurement  of  the  (z,  4>).  The  wafer 
position  itself  determines  r. 

The  Central  Drift  Chamber 

The  CLEO  III  central  drift  chamber  (DR  III)  is  full  of  a  gas  mixture  with  60% 
Helium  and  40%  propane  held  at  about  270  K  and  at  a  pressure  slightly  above  one 
atmosphere.  The  drift  chamber  is  strung  with  array  of  anode  wires  of  gold-plated 
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1  cm 


Figure  2-4:  View  of  the  SVD  III  along  the  beampipe. 


tungsten  of  20  ^m  in  diameter  and  cathode  wires  of  gold-plated  aluminum  tubes 
of  130  fim  in  diameter.  The  anodes  are  kept  at  a  positive  voltage  (about  2000  V), 
and  the  cathodes  are  kept  grounded,  which  provides  an  electric  field  between  the 
anode  and  the  cathode  wires.  Anode  and  cathode  wires  are  often  called  "sense" 
and  "field"  wires  respectively. 

As  a  charged  particle  passes  through  the  DR  III,  it  interacts  electromagneti- 
cally  with  the  gas  molecules  giving  energy  to  the  outer  electrons  which  become  free 
in  a  process  called  ionization.  The  free  electrons  from  the  ionized  gas  molecules 
drift  in  the  electric  field  toward  the  nearest  anode  wire.  As  the  electrons  get  close 
to  the  anode,  the  electric  field  becomes  very  strong  which  causes  an  avalanche 
as  further  ionization  is  induced.  The  result  of  the  avalanche  is  a  large  number  of 
electrons  collapsing  upon  the  sense  wire  in  a  very  short  amount  of  time  (less  than 
one  nanosecond).  When  this  happens  to  a  sense  wire,  we  say  that  there  is  a  "hit". 
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The  current  on  the  anode  wire  from  the  avalanche  is  amplified  and  collected  at 
the  end  of  the  anode  wire.  Both  the  amount  of  charge  and  the  time  it  takes  it  to 
move  to  the  end  of  the  detector  are  measured.  A  calibration  of  the  drift  chamber 
is  used  to  convert  the  amount  of  charge  to  a  specific  ionization  measurement  of  the 
incident  particle.  A  calibrated  drift  chamber  can  also  convert  the  time  to  roughly 
measure  the  position  along  the  sense  wire  where  the  charge  was  deposited. 

The  wires  are  strung  along  the  z  direction.  About  2/3  of  the  outer  part  of 
the  drift  chamber  (the  farthest  part  from  the  interaction  point)  is  strung  in  with  a 
slight  angle  (about  25  miliradians)  with  respect  to  the  z  direction  to  help  with  the 
z  measurement.  Wires  strung  in  the  z  direction  are  called  "axial"  wires,  while  those 
that  are  strung  at  a  slight  angle  are  called  "stereo"  wires. 

The  DR  III  consists  of  an  inner  stepped  section  with  16  axial  layers,  and  an 
outer  part  with  conical  endplates  and  31  small  angle  stereo  layers.  There  are  3 
field  wires  per  sense  wire  and  they  approximately  form  a  1.4  cm  side  square.  The 
drift  resolution  is  around  150  fim  in  r  -  (j)  and  about  6  mm  in  z.  All  wires  are  held 
at  sufficient  tension  to  have  only  a  50  /im  gravitational  sag  at  the  center  (z  —  0). 
There  are  1696  axial  sense  wires  and  8100  stereo  sense  wires,  a  9796  total. 
2.2.3     Crystal  Calorimeter 

The  CLEG  Crystal  Calorimeter  (CC)  is  composed  of  7784  thallium- doped 
Csl  crystals.  Each  crystal  is  30  cm  long  (16.2  radiation  lengths)  with  5cm  x  5cm 
square  front  face.  The  crystals  absorb  any  incoming  electron  or  photon  which 
cascades  into  a  series  of  electromagnetic  showers.  The  electronic  system  composed 
of  4  photo-diodes  present  at  the  back  of  each  crystal  are  calibrated  to  measure  the 
energy  deposited  by  the  incoming  particle.  Other  incoming  particles  other  than 
photons  and  electrons  are  partially,  and  sometimes  fully,  absorbed  by  the  crystal 
giving  an  energy  reading. 
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The  CC  is  arranged  into  a  barrel  section  and  two  endcaps,  together  covering 
95%  of  the  solid  angle.  The  CC  barrel  section  is  unchanged  since  CLEO  II;  the 
endcaps  have  been  rebuilt  for  CLEO  III  to  accommodate  the  new  CESR  interaction 
region  quadrupoles.  The  barrel  detector  consists  of  an  array  of  6144  crystals,  128 
in  (j)  and  48  in  2,  arranged  in  an  almost- projective  barrel  geometry.  That  is,  the 
crystals  are  tilted  in  z  to  point  to  a  few  cm  away  from  the  interaction  point,  and 
there  is  also  a  small  tilt  in  4>.  The  CC  barrel  inner  radius  is  1.02  m,  outer  radius 
is  1.32  m,  and  the  length  in  z  at  the  inner  radius  is  3.26  m.  It  covers  the  polar 
angle  range  from  32  to  148  degrees.  The  barrel  crystals  are  tapered  towards  the 
front  face  (there  are  24  slightly  different  tapered  shapes),  the  endcap  crystals  are 
rectangular,  but  shaved  near  the  outer  radius  to  fit  in  the  container.  The  CC 
endcaps  consist  of  two  identical  end  plugs,  each  containing  820  crystals  of  square 
cross-section,  aligned  parallel  to  the  beam  line  (not  projective).  There  are  60 
crystals  in  the  "fixed"  portion  of  the  "keystone"  piece  of  the  endcap.  760  in  the 
part  that  slides.  The  keystone  is  made  up  of  two  parts,  one  on  top  that  has  12 
crystals  that  for  mechanical  removal  reasons  is  separate  from  a  container  holding  48 
crystals.  The  endcap  extends  from  0.434  m  to  0.958  m  in  r.  The  front  faces  are  z  — 
±1.308  m  from  the  interaction  point.  It  covers  the  polar  angle  region  from  18  to  34 
degrees  in  +z,  and  146  to  162  in  -z. 

The  photon  energy  resolution  in  the  barrel  (endcap)  is  1.5%  (2.5%)  for  5 
GeVphotons,  and  goes  down  to  3.8%  (5.0%)  for  0.1   GeVphotons. 
2.2.4     Ring  Imaging  Cherenkov  Detector 

The  Ring  Imaging  Cherenkov  (RICH)  detector  [23]  is  a  new  detector  subsys- 
tem for  CLEO  III.  It  replaces  the  CLEO  H.V  time  of  flight  system  designed  to 
measure  particles'  velocities. 
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Cherenkov  radiation  occurs  when  a  particle  travels  faster  than  the  speed  of 

light  in  a  certain  medium, 

■ 
v  >  c/n.  (2.1) 

Where  v  is  the  velocity  of  the  particle,  c  is  the  speed  of  light  in  vacuum,  and  n 
is  the  index  of  refraction  of  the  medium  the  particle  is  traveling  in.  The  charged 
particle  polarizes  the  molecules  of  the  medium,  which  then  turn  back  rapidly 
to  their  ground  state,  emitting  radiation.  The  emitted  light  forms  a  coherent 
wavefront  if  v  >  c/n  and  Cherenkov  light  is  emitted  under  a  constant  Cherenkov 
angle,  5,  with  the  particle  trajectory  forming  a  cone  of  light.  The  cone  half-angle  is 
given  by  the  Cherenkov  angle  which  is, 

cos^-^-1-.  (2.2) 

vn      pn 

If  the  radiation  angle,  6,  is  measured,  the  speed  of  the  incident  particle  is  known. 

This  measurement,  combined  with  the  momentum  measurement  from  the  tracking 

system,  gives  a  measurement  of  the  particles  mass,  and  can  be  used  in  particle 

identification. 

The  threshold  velocity  at  which  Cherenkov  radiation  is  emitted  is  vmin  =  j[. 

When  a  particle  traveling  at  the  threshold  velocity  transverses  the  medium  a  very 

small  cone  with  S  »  0  is  produced.  The  maximum  emission  angle  occurs  when 

Vraax  =  c  and  is  given  by 

cos8max  =  -.  (2.3) 

n 

The  RICH  (see  Figure  2-5)  consists  of  30  modules  in  phi,  0.192  m  wide  and 
2.5  m  long.  The  detector  starts  at  a  radius  of  0.80  m  and  extends  to  0.90  m.  Each 
module  has  14  panes  of  solid  crystal  LiF  radiator  at  approximately  0.82  m  radius, 
0.192  m  wide,  0.17  m  long,  1  cm  thick.  Inner  separation  between  radiators  is 
typically  50/mi.  The  LiF  index  of  refraction  is  n  =  1.5.  The  radiators  closest  to  z 
=  0  in  each  module  have  a  45  degree  sawtooth  outer  face,  to  reduce  total  internal 
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reflection  of  the  Cherenkov  light  for  normal  incident  particles  (see  Figure  2-6).  The 
radiators  are  followed  by  a  15.7  cm  (radial)  drift  space  filled  with  pure  N2.  The 
drift  space  is  followed  by  the  photodetector,  a  thin-gap  multiwire  photosensitive 
proportional  chamber. 


Figure  2-5:  The  RICH  detector  subsystem. 


Figure  2-6:  The  two  kinds  of  RICH  LiF  radiators.  For  normal  incidence  particles 
(z  «  0)  a  sawtooth  radiator  is  necessary  to  avoid  internal  reflection. 


With  this  index,  particles  in  the  LiF  radiator  with  beta  =  1  produce 
Cherenkov  cones  of  half-angle  cos_1(l/n)  =  0.84  radians.  With  a  16  cm  drift 


28 

space,  this  produces  a  circle  of  radius  13  cm.  The  RICH  is  capable  of  measuring 
the  Cherenkov  angle  with  a  resolution  of  a  few  miliradians  (see  Figure  2-7).  This 
great  resolution  allows  for  good  separation  between  pions  and  kaons  up  to  about  3 
GeV  as  Figure  2-8  shows. 


3         4         5 
Radiator  row 


Figure  2-7:  Cherenkov  angle  resolutions  per  track  as  a  function  of  radiator  row  for 
Bhabha  events.  Row  1  corresponds  to  the  two  rows  closest  to  z  =  0,  etc. 

2.2.5     Muon  Detectors 

The  muon  detectors  (MU)  are  the  most  external  subsystem  of  the  CLEO  III 
detector.  They  remain  unchanged  from  CLEO  II,  and  are  composed  of  plastic 
proportional  tubes  embedded  in  the  magnet  iron  return  yoke.  They  cover  85% 
of  the  47r  solid  angle  (roughly  30-150  degrees  in  polar  angle).  If  a  series  of  hits  is 
detected  in  the  muon  chamber  layers  they  most  likely  correspond  to  muons  because 
other  particles  are  blocked  by  the  iron.  Besides  detecting  muons,  the  heavy  iron 
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Figure  2-8:  Pion  fake  rate  as  a  function  of  particle  momentum  for  kaon  efficiency  of 
80%  (circles),  85%  (squares)  and  90%  (triangles). 

of  the  return  yoke  protects  the  inner  subsystems  of  the  CLEO  III  detector  from 
cosmic  ray  background  (except  for  cosmic  ray  muons  of  course). 

There  are  three  planes  of  chambers  in  the  barrel  section,  arranged  in  8  octants 
in  (p.  The  plastic  barrel  planes  lie  at  depths  of  36,  72,  and  108  cm  of  iron  (at 
normal  incidence),  corresponding  to  roughly  3,  5,  and  7  hadronic  interaction 
lengths  (16.8  cm  in  iron)  referred  to  as  DPTHMU.  There  is  one  plane  of  chambers 
in  each  of  the  two  endcap  regions,  arranged  in  4  rough  quadrants  in  d>.  They  lie 
at  z  -±  2.7  m,  roughly  covering  the  region  0.80  <  |  cos(0)|  <  0.85.  The  planar 
tracking  chambers  use  plastic  proportional  counters  at  about  2500  V  with  drift  gas 
of  60%  He,  40%  propane,  identical  to  (and  supplied  by  the  same  system  as)  the 
drift  chamber  gas.  Individual  counters  are  5  m  long  and  8.3  m  wide,  with  a  space 
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resolution  (along  the  wire,  using  charge  division)  of  2.4  cm.  The  tracking  chambers 
are  made  of  extruded  plastic,  8cm  wide  by  1  cm  thick  by  5  m  long,  containing  eight 
tubes,  coated  on  3  sides  with  graphite  to  form  a  cathode,  with  50  /zm  silver-plated 
Cu-Be  anode  wires  held  at  2500  V.  The  orthogonal  coordinate  is  provided  by 
8cm  copper  strips  running  perpendicular  to  the  tubes  on  the  side  not  covered  by 
graphite. 
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CHAPTER  3 
ANALYSIS  OF  THE  DATA 

This  work  builds  on  the  techniques  developed  by  previous  CLEO  radiative 

T(1S)  analyses   [18],  [19],  modified  to  be  used  with  CLEO  III  data.  A  new 

technique  based  on  kinematic  fitting  is  developed  to,  together  with  the  new 

RICH  detector,  improve  efficiency  and  particle  identification. 

We  search  for  radiative  T(1S)  decays  in  the  modes  T(1S)  -»  77r+7r~ ^K+K" , 

and  7pp.  The  e+e~  collision  data  has  both  resonant  events,  where  the  e+e~ 

annihilate  to  give  a  T(1S)  ,  and  continuum  events,  where  the  e+e~  collision  does 

not  give  a  T(1S)  .  To  be  sure  we  are  observing  T(1S)  and  not  a  continuum  process, 

the  continuum  must  be  subtracted  by  using  a  pure  source  of  correctly  scaled 

continuum  events.  Pure  continuum  data  can  be  obtained  by  operating  CESR  at  an 

energy  different  form  the  T(1S)  mass.  After  subtracting  the  underlying  continuum, 

we  examine  the  di-hadron  invariant  mass  spectrum  in  search  of  resonances.  We 

determine  the  spin  and  production  helicity  (the  projection  of  the  spin  on  the 

momentum  vector  at  production  time)  of  any  found  resonances  by  examining  the 

photon  and  hadron  angular  distributions. 

3.1     Data  Sample 

The  analysis  presented  here  is  based  on  CLEO  III  data.  Throughout  this 

document,  and  unless  otherwise  stated  we  use  77  luminosity  (the  luminosity  types 

used  by  CLEO  are  defined  at  the  end  of  Section  2.1).  We  prefer  to  use  the  77  over 

the  Bhabha  luminosity  because  the  resonant  process  T(1S)  — >  e+e~  artificially 

increases  the  reported  Bhabha  luminosity  by  about  3%  in  T(1S)  data,  This  extra 

contribution  would  need  to  be  accounted  for  when  doing  a  continuum  subtraction. 

Choosing  the  77  luminosity  avoids  this  complication. 
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The  CLEO  III  data  is  divided  into  numbered  sets.  Sets  18  and  19  have  a 
luminosity  of  1.13  ±  0.02/fe-1  in  the  beam  energy  range  4.727-4.734  GeV.  This 
data,  which  we  call  the  T(1S)  data  (or  simply  the  IS  data),  has  both  resonance 
e+e-  _*  T(1S)  and  continuum  events.  We  take  the  number  of  resonant  events 
from  [24],  JVT(iS)  =  (2-1  ±  0.1)  x  107.  This  number  contrasts  with  the  previous 
generation  measurement  of  CLEO  II,  where  At(is)  «  0.15  x  107  were  available. 

In  datasets  and  18  and  19  there  are  also  0.192  ±  0.004/6"  *  taken  below  the 
T(1S)  beam  energy  (4.714-4.724  GeV).  This  data,  which  we  call  the  T(lS)-off  data 
(or  simply  the  lS-off  data),  has  relatively  low  statistics  and  corresponds  to  purely 
continuum  events. 

To  improve  our  continuum  statistics  we  use  3.49  ±  0.07/6-1  from  datasets 
9,  10,  12.  13,  and  14  of  data  taken  near  the  center-of-mass  energy  of  the  T(4S), 
which  for  our  purposes  is  defined  as  data  with  beam  energy  in  the  5.270-5.300  GeV 
range.  This  set  of  data,  which  we  call  the  T(4S)  data  (or  simply  the  4S  data),  is  a 
source  of  pure  continuum  because  no  T(4S)  —>  BB  resonant  event  can  survive  our 
"cuts"  l    (the  cuts  are  presented  in  Section  3.3). 
3.1.1     Continuum  Subtracted  Distributions 

We  use  the  continuum  data  taken  at  the  T(4S)  energy  to  subtract  the  under- 
lying continuum  present  in  the  T(1S)  data.  This  is  important  because  continuum 
background  processes  like  e+e~   — >  7p  with  p  —>  ir+n',  e+e~  — »  70  with 
0  — »  K+K~,  and  direct  e+e~  — ►  ^h^h"  (we  will  use  the  convention  h  =  7T,  K,  p 
from  now  on),  look  like  the  signal  events  we  are  searching  for.  To  first  order,  the 
cross  section  of  these  continuum  process  scales  like  1/s,  where  s  is  the  square  of 


1  Cuts  are  simply  conditions  that  an  event  must  satisfy  to  be  considered  in  the 
-analysis.  Cuts  are  necessary  to  eliminate  background  that  would  otherwise  make  a 
measurement  difficult  or  even  impossible. 
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the  center  of  mass  energy  of  the  e+e~  system.  Taking  the  luminosity- weighted 
average  beam-energies  of  each  interval,,  and  the  77  luminosities  (see  Table  3-1),  we 
calculate  that  the  T(4S)  data  scaled  down  by  a  factor  of  0.404  2   represents  the  un- 
derlying continuum  in  the  T(1S)  data.  This  is  true  up  to  differences  in  momentum 
distributions  and  phase  space.  The  error  in  the  continuum  scale  factor  is  unknown 
because  the  luminosity  ratio  is  expected  to  have  a  small  but  undetermined  sys- 
tematic error.  We  make  the  somewhat  arbitrary  decision  to  retain  three  significant 
digits  in  the  continuum  scale  factor  because  it  is  sufficient  for  our  purposes  and 
there  are  0.5%  effects  from  second  order  terms  in  the  cross  section  formulae. 

To  eliminate  the  contribution  of  continuum  events  from  a  T(1S)  data  variable 
distribution  (e.g.,  the  invariant  mass  of  two  tracks,  the  photon  angular  distribution) 
we  proceed  as  follows, 

•  Obtain  the  T(4S)  data  distribution  for  the  same  variable. 

•  Efficiency  correct  both  the  T(1S)  and  T(4S)  distributions,  using  a 
GEANT  [25]  based  Monte  Carlo  (MC)  simulation  of  the  detector.  Examples  of 
MC  efficiency  distributions  are  shown  in  Figures  3-3  and   7-3. 


2  If  we  are  to  be  mathematically  strict  we  should  calculate  the  scale  factor  as 

V~^  *~is  run      I    S~~^      *~4S  run 

,c  S1S  run    '    AC  _.,„.  S4S  run 

16  runs  46  runs 

with  obvious  notation.  This  is  equivalent  to  redefining  the  average  energy  as 

J-   \    ^  *-run 


£2      c  a-  m 


runs 


run 


However,  our  energy  intervals  are  sufficiently  narrow  and  this  calculation  does 
not  change  the  last  significant  digit  of  the  original  scale  factor.  Similarly,  taking 
into  account  second  order  terms  in  the  energy  dependence  of  the  cross  section,  like 
the  one  that  appears  as  m2p/s  in  the  explicit  formula  for  the  e+e~    — ♦    jp  cross  sec- 
tion, also  has  an  insignificant  effect  on  the  scale  factor. 
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•  Subtract  the  T(4S)  distribution  from  the  T(1S)  distribution  using  a  T(4S) 
scale  factor  of  0.404. 

We  call  this  set  of  steps  "continuum  subtraction"  by  definition.  For  com- 
pactness, we  call  such  a  distribution  the  "continuum  subtracted  <variable> 
distribution/plot"  or  "<variable>  continuum  subtracted  distribution/plot" .  Except 
for  statistical  fluctuations  and  phase  space  effects,  the  resulting  IS  distribution 
should  not  have  any  contribution  from  continuum  processes  that  scale  as  1/s. 

It  is  important  to  notice  that  any  continuum  subtracted  distribution  is 
efficiency  corrected.  This  means  that  a  fit  to  a  continuum  subtracted  distribution 
(for  example,  the  continuum  subtracted  invariant  mass  distribution)  gives  the 
efficiency  corrected  number  of  events  directly.  Strictly  speaking,  this  number  of 
events  is  only  correct  if  all  the  other  variables  we  cut  on  have  the  same  initial 
distribution  in  data  and  MC,  or  if  the  efficiency  does  not  depend  on  them.  In  this 
note  we  use  "flat  MC" ,  defined  as  MC  that  is  generated  with  a  flat  distribution  in 
the  mass  and  the  helicity  angles  #7,  6*  (these  angles  are  defined  in  the  appendix). 

For  example,  the  number  of  events  obtained  from  the  continuum  subtracted 
invariant  mass  fit  needs  to  be  corrected  to  account  for  the  fact  that  the  helicity 
angle  distributions  are  not  flat  in  data  (see  for  example  Figure  4-9,  and  the 
efficiency  is  highly  dependent  on  these  variables  (see  Figure  7-3).  This  correction  is 
done  in  Section  7.2. 

A  summary  of  the  results  from  this  section  is  shown  in  Table  3-1. 

3.2     Event  Selection 

Event  selection  for  T(15)  — *  ^h+h~  is  straightforward  and  can  be  thought  of 
in  terms  of  three  major  stages. 

First  we  skim  the  data,  keeping  only  those  events  that  contain  exactly  one 
high-energy  photon  and  two  tracks.  Next,  we  require  that  the  total  4-momentum 
of  these  three  elements  be  consistent  with  the  colliding  e+e~  4-momentum.  Finally, 
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Table  3-1:  Summary  of  the  data  used  in  this  analysis.  The  continuum  scale  factor 
is  obtained  using  77  luminosities  because  the  Bhabha  luminosity  is  artificially  high 
diuing  T(1S)  running  due  to  the  process  T(1S)  ->  e+e". 




T(1S) 

T(4S) 

T(lS)-off 

Dataset 

18,  19 

9,  10,  12,  13,  14 

18,  19 

Average  £(,eam(GeV) 

4.730 

5.286 

4.717 

Range  of  ^6eam(GeV) 

4.727-4.734 

5.270  -  5.300 

4.714  -  4.724 

£(e+e-)(/6-1) 

1.20  ±0.02 

3.56  ±  0.07 

0.201  ±  0.004 

JC(77)(/^1) 

1.13  ±0.02 

3.49  ±  0.07 

0.192  ±0.004 

T(1S)  continuum  scale  factor 

1 

0.404 

5.84 

in  the  third  stage,  we  project  the  surviving  data  onto  the  three  different  hadronic 
modes  via  hadron  separation  and  QED  suppression  cuts. 

However,  checking  the  4-momentum  involves  using  the  tracks  masses.  This 
means  that  the  information  from  stage  2  should  somehow  be  useful  in  stage  3.  This 
is  indeed  the  case,  and  the  details  of  how  we  do  it  are  revealed  in  this  section. 
3.2.1     Skim  Cuts 

We  skim  over  the  data  in  the  "hardGam"  subcollection.  The  hardGam 
subcollection  was  developed  with  this  type  of  analysis  in  mind.  For  an  event  that 
passes  the  triggers  3   to  be  classified  as  hardGam  it  must  pass  the  following  cuts, 

•  eGaml  >  0.5 

•  eSh2  <  0.7 

•  eOverPl  <  0.85 

•  eVis  >  0.4 

•  aCosTh  <  0.95  where  "eGaml"  is  the  highest  isolated  shower  energy 
relative  to  the  beam  energy,  "eSh2"  is  the  energy  of  the  second  highest  shower 
relative  to  the  beam  energy,  "eOverPl"  is  the  matched  shower  energy  relative  to 


3  Triggers  are  basic  criteria  that  an  event  must  satisfy  to  recorded  during  the 
data  collection  processes.  Triggers  are  designed  to  get  rid  of  trash  and  noise  and 
reduce  the  size  of  the  data  sample  while  keeping  all  of  the  important  information. 


36 

the  momentum  of  the  track  with  highest  momentum,  "eVis"  is  the  total  energy 
detected  (charged  tracks  are  assumed  to  be  pions)  relative  to  the  center  of  mass 
energy,  and  "aCosTh"  is  the  z  component  of  the  unitary  total  momentum  vector. 

Monte-Carlo  predicts  that  about  75%  of  the  generated  T(1S)  -»  7/i+/r  signal 
passes  the  hardware  and  software  triggers  and  gets  classified  as  hardGam. 

As  mentioned  above,  we  use  this  data  to  make  our  skim.  To  write  an  event 
from  the  hardGam  subcollection  into  our  skim  we  require  the  following  topological 
cuts, 

•  There  are  exactly  two  "good  tracks" :  there  can  be  any  number  of  tracks  that 
are  not  "good  tracks"  but  these  are  not  used  in  the  analysis.  We  define  a  "good 
track"  as  a  track  that  satisfies  the  following  cuts;  drift  chamber  track  ionization 
energy  loss  (dE/dX)  information  is  available,  the  ratio  of  number  of  wire  hits  to 
those  expected  is  between  0.5  and  1.5,  the  pion  fit  has  x2/d.o.f.  <  20  (here  d.o.f. 
stands  for  degrees  of  freedom),  and  the  distance  of  closest  approach  to  the  beam 
spot  in  the  x-y  plane  (called  DBCD)  is  less  than  5  -  3.8P  (mm)  if  P  <  1  GeV/c 
(where  P  is  the  tracks  momentum  in  GeV/c)  and  less  than  1.2  mm  for  tracks  with 
P  >  1  GeV/c.  This  DBCD  cut  is  common  in  the  more  sophisticated  CLEO  II/II.5 
analyses.  It  performs  better  than  a  simple  DBCD  <  5  mm  cut,  because  it  takes 
into  account  the  fact  that  traks  with  higher  momentum  have  a  better  measurement 
of  DBCD  since  they  scatter  less. 

•  There  is  exactly  one  "good  shower" ,  there  can  be  any  number  of  showers 
that  are  not  "good  showers"  but  these  are  not  used  in  the  analysis.  We  define  a 
"good  shower"  as  an  unmatched  shower  with  energy  >  4  GeV. 

These  topological  cuts  are  about  85%  efficient  for  generated  signal  events  that 
have  passed  the  triggers  and  have  been  classified  as  hardGam. 

The  overall  skim  efficiency  is  between  60-65%,  depending  on  the  mode  (see 
Table  3-4). 
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3.2.2     Analysis  Cuts 

After  our  skim  we  call  any  cuts  we  make  "analysis  cuts" .  These  cuts  are  done 
at  analysis  time  and  are  mode  dependent.  As  a  convention,  and  unless  otherwise 
stated,  efficiencies  for  individual  analysis  cuts  are  reported  relative  to  the  events  in 
the  skim  (not  relative  to  the  events  generated). 
4-momentum  Cut 

All  fully  reconstructed  events  should  have  the  4-momentum  of  the  eTe" 
system.  This  constraint  is  usually  implemented  with  a  simple  two-dimensional 
AE-p  box  cut,  where  AE  is  the  difference  between  the  reconstructed  energy  for 
the  event  and  the  colliding  <-:*  e"  energy  (Ecm),  and  p  is  the  magnitude  of  the 
reconstructed  total  momentum  for  the  event.  Typical  values  for  these  cuts  are 
-0.03  <  AE/Ecm  <  0.02  and  p  <  150  MeV/c  (taken  from  [19]). 

The  traditional  AE-p  box  cut  is  somewhat  useful.  However,  it  does  not  take 
into  account  the  correlation  between  the  measured  energy  and  momentum.  Indeed, 
the  signal  lies  in  diagonal  bands  in  the  AE-p  plane,  making  a  box -shaped  cut  not 
optimal  (see  Figures  3-1  a  and   3 -lb). 

We  use  an  alternative  approach  to  the  4-momentum  cut.  After  a.  simple 
substitution,  J57  =  p1  =  ECm  ~  Eh+h-  (where  £7  is  the  photon's  energy,  p7  is 
the  magnitude  of  the  photons  momentum,  Eqm  is  the  energy  of  the  e+e~  system, 
and  Eh+h-  is  the  energy  of  the  hadron  pair),  we  can  write  the  E-p  conservation 
equations  as: 

Ph+h-  +  (Ecm  -  Eh+k-  )?y  ■=  Pcm  (3-1) 

where,  Ph+h-  is  the  di-hadron  momentum,  jo7  is  the  photon's  momentum  unit 
vector,  and  pcu  is  the  momentum  of  the  e+e~  system  (which  is  a  few  MeV  because 
of  the  crossing  angle).  Equation  3.1  is  a  3-constraint  subset  of  the  4-momentum 
constraint  and  has  the  convenient  property  of  avoiding  the  use  of  the  measured 
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photon's  energy,  which  has  non-Gaussian  asymmetric  errors.  It  is  important  to 
notice  that  Equation  3.1  contains  the  di-hadron  energy,  therefore,  it  can  help 
discriminate  between  the  various  particle  hypotheses. 

We  proceed  as  follows.  After  vertexing  the  hadron  pair  using  the  beam  spot 
with  its  error  matrix,  we  calculate  the  photon's  direction  from  the  hadron  pair's 
vertex  and  the  shower  position.  We  then  fit  the  event  to  the  3  constraints  expressed 
in  Equation  (1)  using  the  techniques  outlined  in  [26]  and  cut  on  the  x2  of  the 
3-constraint  fit,  xl-pCO  <  10°-  To  complete  the  4-momentum  requirement,  we 
calculate  AE(h)  =  Eh+h-  +  E1  -  ECm,  where  Eh+h-  is  the  updated  di-hadron 
energy  after  the  constraint,  and  £7  is  the  measured  photon's  energy,  and  require 
-0.050  *  ECm  <  &E(h)  <  0.025  *  ECm-  Furthermore,  we  now  have  available 
X%-P{h)  differences  between  different  particle  hypotheses,  which  help  in  particle 
identification  (ID).  This  is  discussed  in  more  detail  in  the  next  section.  Figure  3-1 
compares  the  performance  of  the  old  and  new  approaches  to  the  E-p  cut. 

At  this  point  it  is  a  good  idea  to  check  that  the  4-momentum  cut  rejects 
background  events  that  make  it  through  our  skim  cuts  (see  Section  3.1).  These 
events  typically  have  one  high  energy  shower,  two  tracks,  and  (an)  additional 
element(s).  One  such  background  is  T(1S)  -»  ^tt'tt0.  Out  of  25000  T(15)  -»  7»?, 
with  7/  -*  7T+7r~7r°  MC  events  4   only  4  survive  our  4-momentum  cut.  We  conclude 
that  our  4-momentum  cut  is  good  at  rejecting  background  events  that  pass  the 
skim  cuts  but  have  additional  elements  such  as  an  extra  photon,  7T°,  pair  of  tracks, 
etc. 

Hadron  Separation 

We  define  Ax/D(/ij  -  h2)  for  the  particle  hypotheses  hi  and  h2  of  our  charged 
track  pair  (e.g.  Ax/£>(7r  —  K))  as  follows, 


4  Thanks  to  Vijay  Potlia  for  generating  these  events 
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Figure  3-1:  Distributions  for  different  4- momentum  cuts.  Signal  MC  events 
(T(15)  — »  77r+7r~)  are  represented  by  the  black  dots,  and  "background"  MC  events 
(e+e~   -*  -yep,  <j>  -*  K+K~)  appear  as  (red)  triangles.  Plot  a)  has  no  cuts.  Plot  b) 
has  the  old  4- momentum  box-cut.  Plot  c)  has  the  new  4- momentum  cuts.  Plot  d) 
has  the  new  4- momentum  cuts,  and  also  a  cut  defined  as  xl-p(7r)  ~  Xe-p(^)   <  0. 
The  particle  ID  potential  of  the  newly  available  XE-P(h)  is  evident. 


A.X/D^l    -  h)  =      °dE/dx(ht)    -  VdE/dxiH)  +  °iE/dx(hl)  -  adE/dx(h2) 

-2\og(CRjcH{hi))  +  2log(£mCH(K))  (3-2) 

-2\og{CmcH(fh))  +21og(£  mcniK)). 

where  the  idea  is  to  combine  the  dE/dX  and  RICH  information  into  one  number. 
Pairs  of  tracks  with  AxjZ)(/i1  —  /i2)  <  0  are  more  likely  to  be  of  type  hi  than  of  type 

h2. 

In  practice,  we  only  add  a  tracks  RICH  information  if  its  momentum  is  above 
the  Cherenkov  radiation  threshold  for  both  mass  hypotheses  and  there  are  at  least 
3  photons  within  3a  of  the  Cherenkov  angle  for  at  least  one  of  the  mass  hypothesis. 
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We  also  require  that  both  hypotheses  were  actually  analyzed  by  RICH  during 
pass2. 

In  addition  to  RICH  and  dE/dX,  and  as  hinted  in  the  previous  section,  the 
difference  in  xi-pM  fr°m  tne  constraint  expressed  in  Equation  3.1  can  help  the 
particle  ID  (see  Figures  3-lc  and  3- Id).  We  define, 

AxL#i  -  fca)  =  xl-pCn)  "  xl-pW-  (3-3) 

Events  with  Ax|_p(/ii  -  h2)  <0  are  more  likely  to  be  of  type  hi  than  of  type  h2. 

In  this  analysis,  to  select  hi  and  reject  h2  the  default  cut  is  simply  Ax/£>(hi  - 
h2)  <  0.  This  simple  cut  is  highly  efficient,  has  low  fake  rates,  and  is  sometimes 
sufficient.  However,  out  of  the  six  possible  cases  when  one  hadron  fakes  another, 
there  are  three  important  cases  where  it  pays  off  to  also  use  Ax|_p(hi  -  h2) 
together  with  Ax'fD(/ii  -  h2)  in  an  optimal  way; 

1.  -it  background  to  K.  This  background  comes  from  the  continuum  process 
e+e_  —*  7p,  p  —*  7r+'/r~. 

2.  7r  background  to  p.  Again,  this  background  comes  from  the  continuum 
process  e+e~  — >  7/9,  p  — *  7:^7r~. 

3.  K  background  to  p.  This  background  comes  from  the  continuum  process 
e+e~  — ►  7^,  <f>  —*  K+K~. 

In  other  words,  the  important  cases  occur  when  the  lighter  mass  hypothesis 
fakes  the  heavier  mass  hypothesis. 

Mathematically,  one  would  expect  that  simply  adding  both  Ax2ID  and  Ax|_p 
together  (like  we  just  did  when  combining  RICH  and  dE/dX),  and  cutting  on  the 
grand  Ax2  is  the  way  to  go.  Unfortunately,  this  simple  approach  fails  because  of 
large  non-mathematical  tails  in  the  individual  xl-P  distributions. 

Instead,  for  each  of  these  three  cases  we  define  the  best  cut  values  (ci,  c2)  in 
AxjD  <  C\  and  Ax|;_p  <  c2  as  those  that  maximize 
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F{ci,cs)  *=• 


IL 


(3.4) 


where  #e  (i?/)  is  the  efficiency  (fake  rate)  of  the  particle  ID  cuts  and  W  is  the 
rough  ratio  of  the  background  to  signal  in  the  data  sample  for  each  case.  Each  W 
can't  be  known  a  priori,  but  a  rough  idea  of  its  value  can  be  obtained  by  doing 
a  first  iteration  of  the  analysis  with,  for  example,  c\  =■  oi  —  0.  We  use  W  —  20, 
W  —  60,  and  W  —  30  for  cases  1-3  respectively.  Figure  3-2  shows  F(ci,  c2)  for  each 
case.  Table  3-2  shows  the  optimized  cut  values  and  their  effect  on  particle  ID. 


Figure  3-2:  Contour  plots  of  F(ci,  c2)  as  defined  in  Equation  3.4  for  different  sig- 
nals and  backgrounds  from  flat  MC.  Top  has  a  K  signal  and  a  n  background.  Cen- 
ter has  a  p  signal  and  a  vr  background.  Bottom  has  a  p  signal  and  a  K  background. 


Quantum  Electrodynamic  Suppression 

Quantum  Electrodynamic  (QED)  background  in  this  analysis  comes  from  the 
abundant  processes  e+e~  — >  7e+e'~  and  e+e"  — *■  7//f//~. 
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Table  3-2:  Cut  values,  efficiencies,  and  fake  rates  for  Axw  <  ci  anc^  &Xe-p  <  c2  in 
flat  MC.  Trie  different  ci,  c2  values  are  chosen  so  F  from  Equation  3.4  is  maximun. 
Efficiency  and  fake  rates  of  each  cut  are  reported  relative  to  events  in  the  skim. 
Statistical  errors  in  the  efficiencies  are  0.1%  or  less.  Errors  in  the  fake  rates  are 
statistical!  and  are  shown  for  completeness  only. 


7r  faking  K 
k  faking  p 
K  faking  p 


Cut  Value 


3 
-1 

-1 


-2 
-2 
3 


Signal  efficiency  (%) 
AxL>    AxL,    Both 


94.2 
99.0 

98.8 


95.2 
99.8 
99.3 


89.9 
98.9 
98.1 


&Xw 


Fake  rate  (%) 

Ze-p 


Ax] 


Both 


1.73  ±0.09  14.5  ±0.3  0.31  ±  0.04 
0.89  ±0.06  1.03  ±0.07  0.03  ±0.01 
1.89  ±0.09     1.09  ±0.07    0.10  ±0.02 


To  reject  e+e    —*  ^e+e    we  require  both  tracks  to  have  a  matched  shower  with 
energy  E  such  that  \E/p(ir)  -  0.95|  >  0.1,  and  to  have  AxjD(h  -  e)  <  0. 

To  reject  efe"  -»  yp,+lT  we  simply  require  DPTHMU  <  5  (DPTHMU 
was  defined  in  Section  2.2.5)  for  both  tracks  in  the  K  and  p  modes  because  particle 
ID  cuts  make  the  pion  (and  therefore  the  muon)  fake  rate  small  (see  Tables  3-2 
and  3-6).  For  the  7r  mode  we  cannot  use  particle  ID  in  a  practical  way  because 
muons  and  pions  have  similar  masses.  Instead,  to  separate  pions  and  muons  we 
use  a  much  stronger  cut  requiring  that  both  tracks  be  within  the  barrel  part  of 
the  muon  chamber  (|  cos(0)|  <  0.7),  both  have  P  >  1  GeV/c  and  both  have 
DPTHMU  <  55  . 

To  improve  the  overall  muon  suppression  cut  efficiency  with  virtually  no 
increase  in  muon  fakes,  we  flag  an  event  as  "not  muonic"  if  any  of  the  tracks 
deposit  more  than  600  MeV  in  the  CC,  This  increases  the  cut  efficiency  by  about 
90S^  in  the  tt  mode  and  makes  the  detector  more  hermetic. 


5  Other  analyses  (for  example  [18])  typically  use  DPTHMU  <  3,  our  CLEO  III 
MC  has  too  many  pion  tracks  with  3  <  DPTHMU  <  5  and  two  few  tracks  with 
DPTHMU  =  0. 


s      ■ 
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3.2.3     Cut  Summary,  Efficiencies,  and  Fake  Rates 

Table  3-3  summarizes  the  cuts  used  in  this  analysis.  Figure  3-3  shows  the 
overall  Monte  Carlo  efficiency  after  all  cuts.  Figure  3-4  shows  the  fake  rates 
according  to  the  MC  and  the  data  for  different  particle  ID  cuts.  The  data  fake 
rates  and  limit  fake  rates,  which  are  measured  for  pions  and  kaons  faking  other 
hadrons,  are  calculated  from  the  p  and  <f>  peaks  in  the  continuum  (see  Figures  4-2 
and   5-1).  Tables  3-4  through  3-6  summarize  the  results  of  Figures  3-3  and  3-4. 


44 


cc 

>> 

33 
a 

«3 
'i. 

X 


r 


CO 

I 


55 

1 

to 

o 

(h 

X 

CD 

CO 

s~— ^ 

00 

ft. 

bO 

T3 

<D 

s-i 

hi 

<m 

O 

-tJ 

JS 

* 

tu 

1 

o 

.23 

te 

i 

C 

II 

o 

o 

CC 

■C 

-ts 

>. 

9) 

c 

"3 

i- 

ce 

o 

0) 

ts 

^ 

i—H 

c 

"S  § 

CO 

1 

cp   ce 

^ 

.« 

J  u 

y 

+ 

eg     H 
s-i     cS 

■Si 

2?  .-C 

T3 

T 

CO 

o 
o 

QO 

CO 

C 
u 

CO 
•  1— < 
CO 

5^> 

CO 

I— t 

-t-2 
C 

u 

a 

co     O 

CO 

o 

12 

a.o 

fu 

ci 

c 

CO 

£  A 

C 
< 

'H 

+3     C? 

c   a 
w3 

o5  J3 

C 

-c  .-c 

G 

O 

•  t— i 

-fc? 

•— 1 

CO 

j    ,| 

c 
o 

'3 

c6 

T» 

o 

"5b 

$ 

05 

o 

> 

ce 

o 

0 

-4-3 

cS 

Q 

45 


Table  3-4:  Efficiencies  in  %  for  cuts  (as  outlined  in  Table  3-3)  for  fiat  signal  MC. 
Efficiencies  in  the  second  group  are  reported  relative  to  the  number  of  candidates 
that  make  it  to  the  skim.  The  third  part  of  the  table  shows  the  overall  reconstruc- 
tion efficiency.  Statistical!  errors  are  0.1%  or  less. 


Skim  cuts 

7T 

K 

P 

Hardware  Trigger 

Software  Trigger 

hardGam 

Topological  (ntracks  =  2  &  £7  >  AGeV) 

89.7 
96.5 
78.4 
73.3 

89.0 
96.2 
79.9 
67.6 

90.9 
96.6 
76.5 
72.3 

Overall  Skim  efficiency 

64.3 

60.3 

62.4 

Analysis  cuts 

7T 

K 

P 

4-momentum 

QED  e+e"  — »  /ye+e"  suppresion 

QED  fi+H~  —*  7At+/-i~  suppresion 

Hadron  separation 

98.6 
93.9 
74.7 
97.1 

98.5 
87.4 
93.0 
89.0 

99.0 
93.1 
98.3 

97.8 

Overall  analysis  efficiency 

66.9 

79.1 

89.1 

Overall  reconstruction  efficiency 

43.0 

47.6 

55.6 

Table  3-5:  Efficiencies  in  %  for  cuts  (as  outlined  in  Table  3-3)  for  flat  4S  MC.  Effi- 
ciencies in  the  second  group  are  reported  relative  to  the  number  of  candidates  that 
make  it  to  the  skim.  The  third  part  of  the  table  shows  the  overall  reconstruction 
efficiency.  Statistical!  errors  are  0.1%  or  less. 


Skim  cuts 

7T 

K 

P 

Hardware  Trigger 

Software  Trigger 

hardGam 

Topological  (ntracka  =  2  &;  E1  >  AGeV) 

88.4 
94.6 
79.2 
75.2 

88.0 
94.6 
80.5 
69.5 

89.0 
94.8 
77.9 
75.8 

Overall  Skim  efficiency 

66.9 

62.7 

66.6 

Analysis  cuts 

7T 

K 

P 

4-momentum 

QED  e+e~  — ►  ^e+e~  suppresion 

QED  ix+n~  — *  7/if  fjT  suppresion 

Hadron  separation 

97.1 
93.4 
75.8 
98.0 

97.0 
89.1 
92.5 
91.4 

97.4 

.93.2 

98.1 

96.6 

Overall  analysis  efficiency 

67.4 

80.2 

87.3 

Overall  reconstruction  efficiency        |  45.0 

50.3 

58.2 
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Table  3-6:  Final  efficiencies  and  fake  rates  after  all  cuts  in  %.  Statistical  errors 
for  efficiencies  are  0.1%  or  less.  Statistical  errors  for  fake  MC  rates  are  shown  for 
completeness  only.  The  IS  DATA  corresponds  to  the  IS  off  resonance  data  sample. 
The  MC  was  generated  flat. 


7T  CUtS 

K  cuts 

p  cuts 

IS*  MC 

43.0 

0.14  ±0.01 

<  0.007 

4S  7T  MC 

45.0 

0.29  ±  0.02 

<0.01 

IS  p  DATA 

... 

<  1.2 

<0.2 

4S  p  DATA 

-    ■ 

0.20  ±  0.06 

<0.06 

1SXMC 

1.27  ±0.03 

47.6 

<  0.02 

4S  K  MC 

0.92  ±  0.04 

50.3 

<0.06 

IS  <j>  DATA 

<3.8 

- 

<2.0 

4S  4>  DATA 

4.14  ±0.69 

- 

<  0.4 

ISpMC 

0.08  ±  0.02 

0.34  ±0.04 

55.6 

4S  p  MC 

0.05  ±  0.01 

0.37  ±0.01 

58.2 
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Figure  3-3:  Final  efficiencies  for  each  mode  as  a  function  of  invariant  mass  for  the 
IS  (solid)  and  4S  (dashed)  Monte  Carlo  data.  The  MC  was  generated  with  a  flat 
angular  distribution. 
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Figure  3-4:  Fake  rates  for  IS  (hollow  circles)  and  4S  (hollow  squares)  according  to 
flat  MC.  The  (red)  downward  pointing  triangle  is  obtained  using  the  p  peak  in  45 
data. 


CHAPTER  4 
EXCLUSIVE  RADIATIVE  DECAY  T(15)  -  Tt+iT 

In  Figure  4-2  the  k+tt"  invariant  mass  plot  is  shown  for  both  IS  and  4S  data. 
Figure  4-3  shows  the  continuum  subtracted  ir+Tr~'  invariant  mass  distribution  (as 
defined  in  Section  3.2.1)  with  the  most  likely  statistical  fit  overlayed  (which  is 
described  in  the  next  section).  The  number  of  events  within  ir  of  the  p  region 
(0.62  -  0.92  GeV/c2)  left  after  the  continuum  subtraction  is  200  ±  300,  and  50  of 
these  belong  to  the  /o(980)  low- mass  tail. 

4.1     Robustness  of  The  Mass  Distribution 

In  [27]  the  decay  T(1S)  -»  77r°7r°  is  analyzed,  and  it  is  shown  how  the  analysis 
stream  warps  the  shape  of  the  reconstructed  resonance.  This  effect,  which  arises 
because  of  the  particular  n°  behavior,  raises  problems  when  fitting  the  invariant 
mass  distribution. 

In  Section  3.2.1  we  claimed  that  if  the  data  and  MC  had  the  same  91  and 
#7r  distributions,  the  fit  to  the  continuum  subtracted  invariant  mass  distribution 
automatically  gives  the  correct  efficiency  corrected  number  of  events. 

Here  we  test  this  claim.  To  this  end,  we  generate  10000  T(1S)  -»  7/2(1270), 
with  /2(1270)  -+  7r+7T  with  flat  91  and  6V  distributions.  We  treat  this  MC 
as  data  and  carry  out  the  first  two  steps  of  the  continuum  subtraction  process. 
The  resulting  peak  has  a  mass  of  1.278  ±  0.002  GeV/c2  and  a  width  of  0.193  ± 
0.006  GeV/c2,  consistent  with  the  generated  mass  and  width  of  1.275  GeV/c2  and 
0.185  GeV/c2.  More  importantly,  the  number  of  reconstructed  events  from  the  fit  is 
10040  ±  180,  which  is  consistent  with  the  number  of  generated  events. 

Figure  4-1  shows  the  reconstructed  and  efficiency  corrected  events,  a  fit  to 
them,  and  the  generated  events. 
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We  conclude  that  there  is  no  warping  of  the  mass  distribution,  and  that  the 
analysis  stream  behaves  like  we  expect  when  obtaining  the  efficiency  corrected 
number  of  events  from  data. 

■    Reconstructed  and  Eff.  Corrected 
r— 1  MC  Generated 
2400 |    i   i    i    i    |    i    i    i    i    |   r i   i    i    |    i 


T  '   '    '    '    I" 


1 


I    i    i    i    i   I 


1.00  1.10  150  1.30  1.40  1.50 

7t+tt"*  Invariant  Mass 


Figure  4-1:  Reconstructed  events  and  efficiency  corrected  events,  a  fit  to  them,  and 
the  generated  events  for  T(1S)  -*  7.M1270),  with  /2(1270)  ->  tt+tt"  with  flat  07 
and  6-V  distributions. 


4.2     Statistical  Fit  of  the  Invariant  Mass  Distribution:  Signal  Areas  and 

Their  Significance 

The  results  of  this  section  are  summarized  in  Table  4-1. 

Figure  4-3  shows  possible  signals  for  for  T(15)  — >  7/o(980),  T(15)  --> 
7/2(1270)  and  T(15)  -»  7/4(2050).  To  determine  the  number  of  events  in  each 
signal  we  fit  the  invariant  7r+7r~  mass  continuum  subtracted  distribution  with  three 
spin-dependent,  relativistic  Breit-Wigner  line  shapes.  The  masses  and  widths  are 
allowed  to  float,  except  for  the  width  of  the  /4(2050)  which  has  a  very  large  error 
if  allowed  to  float  and  is  set  to  its  PDG  value.  The  PDG  values[16]  for  the  mass 
and  width  of  the  three  resonances  are  m/0(980)  =  980  ±  10  MeV/c2,  T^ggo)  = 
70  ±  15  MeV/c2,  m/2(1270)  =  1275.4  ±  1.2  MeV/c2,  r/2(1270)  =  185.1^  MeV/c2, 
m/4(2050)  =  2034  ±  11  MeV/c2,  and  r/4(205o)  =  222  ±  19  MeV/c2. 
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To  measure  the  statistical  significance  of  each  signal  we  do  multiple  x2  fits 
fixing  the  signal  area  to  different  values  while  letting  the  mass  and  width  of 
the  signal  whose  significance  is  being  measured  to  float  within  2a  of  the  PDG 
values.  At  the  same  time  the  rest  of  the  fit  parameters  are  fixed  to  the  values  that 
originally  minimized  the  x2  ■  We  assign  each  of  these  multiple  fits  a  probability 
proportional  to  e~x  I'1  and  then  normalize.  We  calculate  the  chance  of  the  signal 
being  due  to  a  random  fluctuation  by  adding  the  normalized  probabilities  for  the 
fits  with  a  negative  or  0  signal.  This  method  fails  for  the  highly  significant  /2(1270) 
signal  because  the  e~x  /2  value  of  fits  with  negative  or  0  /2(1270)  signal  is  of  the 
order  of  10-54  and  our  computing  software  can  only  handle  numbers  as  small  as 
10-45.  For  completeness  we  state  that  the  significance  of  this  signal  is  <  10-45. 

To  measure  the  upper  limit  for  the  fj (2220)  we  also  do  multiple  x2  fits 
for  different  fixed  signal  values,  while  keeping  its  mass  and  width  constant  at 
mfj{2220)  =  2.234  GeV/c2  and  r/j(222o)  '=  17  MeV/c2  as  in  [17].  The  resulting 
probability  plot  is  shown  in  Figure  4-5. 

We  find  clear  evidence  for  the  /2(1270),  evidence  for  the  /o(980)  and  weak 
evidence  for  the  /4(2050).  We  also  put  a  90%  confidence  level  upper  limit  on 
/j(2220)  production.  Fit  results  are  shown  in  Table  4-1. 

Table  4-1:  Results  for  T(15)    -»    77r+7T.  The  branching  fractions  of  /2(1270)  and 
/4(2050)  are  taken  from  the  PDG  [16].  Errors  shown  are  satistical  only. 


Mode 

Area 

B.F.  (10-5) 

Significance 

7/o(980),/o(980)^7r+7r- 

340f™~ 

1  fi+U7 

8.3  x  10-6  (4.3a) 

7/2(1270) 

1230  ±  100 

10.4  ±0.8 

<  10-45    (>  14a) 

7/4(2050) 

85  ±  30 

3.6  ±1.3 

5.2  x  10-3  (2.6a) 

7/j(2220),  /j(2220) ->  tt+tt- 

<  13 

<  6.2  x  10-2 

- 

4.3     Angular  Distribution  of  The  Signal 

The  helicity  angle  distributions  of  9n  and  #7  are  defined  and  described  in  the 
appendix.  In  this  section  we  first  obtain  the  helicity  angular  distributions  of  the 
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/o(980),  /2(980),  and  /4(2050)  and  then  fit  them  to  the  predictions  of  the  helicity 
formalism. 

In  practical  terms,  obtaining  the  helicity  angle  distribution  of  a  particular 
resonance  from  data  consists  of  two  steps.  First,  we  choose  an  invariant  mass 
interval  around  the  resonance  peak  to  select  events  from  the  resonance  and  obtain 
a  helicity  angular  distribution  which  has  both  signal  and  background  events. 
Second,  we  subtract  the  contribution  to  the  helicity  angular  distribution  of  the 
background  events  in  the  chosen  mass  interval  to  obtain  what  we  want;  the  helicity 
angle  distribution  of  the  resonance. 

Choosing  the  mass  interval  is  not  a  trivial  thing.  If  its  too  wide  there  will 
be  too  much  background,  and  if  its  to  narrow  there  will  not  be  enough  signal. 
To  choose  the  optimum  mass  interval  we  need  to  know  how  much  signal  and 
background  we  are  selecting.  Therefore,  the  two  steps  described  in  the  previous 
paragraph  are  related.  How  we  deal  with  this  is  revealed  in  the  next  two  sections. 
4.3.1     Optimum  Mass  Interval 

As  described  above,  the  first  step  in  obtaining  the  helicity  angle  distribution 
for  a  resonance  is  choosing  an  invariant  mass  interval  to  select  events  from  such  a 
resonance.  A  standard  IT  (which  corresponds  to  1.6cr  for  a  spin  0  Breit-Wigner) 
cut  around  the  mean  mass  of  the  resonance  can  be  chosen  as  a  "standard"  interval. 
We  could  proceed  this  way,  but  in  our  case  because  of  the  large  subtractions 
involved  when  obtaining  the  angular  distribution,  a  considerable  increase  statistical 
significance  of  each  bin  in  the  helicity  angular  distribution  can  be  achieved  by 
choosing  the  mass  interval  carefully  (see  the  last  column  of  Table  4-2). 

Let's  consider  for  example  the  /o(980)  angular  distribution.  We  begin  with 
the  T(IS')  and  the  T(AS)  efficiency  corrected  distributions.  Before  the  continuum 
subtraction  each  bin  in  the  /o(980)  angular  distribution  has  contributions  from 
the  high  end  mass  tail  of  the  p  and  the  low  mass  tail  of  the  /2(1270).  In  order  to 
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- 

get  the  final  angular  distribution,  both  of  these  contributions  are  taken  away  by 
first  doing  a  continuum  subtraction  using  the  scale  factor  a  =  0.404,  and  then 
by  subtracting  the  /2(1270)  distribution  outside  the  mass  interval  scaled  by  an 
appropriate  scale  factor  (5  equal  to  the  ratio  of  /2(1270)  inside  the  mass  interval 
and  outside  of  it  1  .  After  all  this,  the  contributions  to  the  bins  in  the  final  angular 
distribution  are, 

N/o(980)  =  ?Vo(980)  +  np(T(lS))  +  ™/2(1270)  ~  anp(T(4S))  ~  ^/2(1270)  ~  nM980)         (^l) 

where  nMm{nMim))  is  the  number  of  /0 (980) (/2(  1270))  obtained  by  integrating 
the  fitted  spin-dependent,  relativistic  Breit-Wigner  function  inside  the  mass 
interval,  np{T(iS))(np{y(4S)))  is  the  number  of  p's  inside  the  mass  interval  from 
continuum  events  at  the  T(1S)(T(4S))  energy,  and  n'/?(1270)  is  the  number  of 
/2(1270)  outside  the  mass  interval  being  used  to  subtract  the  contribution  of 
/2(1270)  to  the  /O(980)  angular  distribution.  Following  this  last  definition,  /3  = 

«A(i2TO)/n/a(1270)' 

Each  of  these  terms  has  an  associated  error.  Assuming  that  An  =  y/n,  that 

the  efficiency  correction  has  infinite  statistics,  and  ignoring  the  errors  on  the 
continuum  scale  factor  a,  the  overall  error  is, 


AiV/o(980)  =  ^0(980)  +  ("  +  a2K(T(4S))  +  (2  +  P  +  P2)nf2(1270) 

We  arrive  to  the  conclusion  that  the  mass  interval  (m  -  Am,  fh  +  Am)  which 
produces  the  helicity  angle  distribution  with  smallest  relative  bin-errors  is  the  one 


1  Strictly  speaking  this  is  correct  only  to  first  order.  The  /2(1270)  distribution 
itself  has  a  small  contamination  from  /0(980)  and  /4(2050).  This  effect,  which  we 
call  cross-contamination,  is  ignored  in  Equation  4.1.  Later,  in  section  4.2.2,  we  will 
show  how  to  eliminate  this  small  cross-contamination  using  all  the  resonances.  The 
difference  in  how  we  actually  get  the  helicity  distributions  and  Equation  4.1  has  an 
insignificant  effect  when  calculating  the  optimal  mass  interval 


> 
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that  maximizes, 


FM98o){fn,  Am)  = 


n/o(980) 


(4.2) 


V^/o(980)  +  5  +  a2)np{r(4S))  +  (2  +  P  +  /32)n/2(mo) 
The  plot  of  F/o(98o)(m,  Am)  is  shown  in  figure  4-6. 

This  same  technique  can  be  applied  to  the  /2(1270)  and  /4(2050)  resonances. 
Results  are  shown  in  Table  4-2. 

Table  4-2:  Mean  masses,  widths  in  GeV/c2  and  inverse  of  the  average  bin  relative 
error  (F)  from  background  subtractions  for  the  angular  distribution  of  different 
resonances.  Standard  mean  masses  and  widths,  corresponding  to  ir,  are  taken 
from  the  fit  in  Figure  4-3,  and  are  labeled  with  the  subscript  "s",  while  those  that 
maximize  F  are  labeled  with  the  subscript  "m" .  The  last  column  shows  the  factor 
by  which  the  effective  statistics  increase. 


Resonance 

(ms.  Ams) 

F(m„Am,) 

(mm,  Amm) 

F(mm,Amm) 

F2(mm,&rnm) 
F2(ms,Ams) 

/o(980) 

/2(1270) 

/4(2050) 

(0.970,6.070) 
(1.270,0.120) 
(2.120, 0.220) 

5.7 

18.3. 

3.2 

(0.985, 0.060) 
(1.590,0.420) 
(2.240,0.250) 

6.0 

20.6 

3.4 

1.11 

1.27 
1.13 

4.3.2     Background  Subtraction 

The  continuum  subtraction  of  the  helicity  angle  distribution  is  defined  in 
Section  2.1.  The  subtraction  of  the  tails  from  other  resonances  requires  a  closer 
look.  Let  us  call  the  continuum  subtracted  helicity  angle  distribution  of  the  events 
in  the  /o(980),  /2(1270),  and  /4(2050)  mass  intervals  c/o(98o),  C/2(i27o),  and  c/4(205o) 
respectively.  Let  us  call  the  helicity  angle  distribution  of  the  events  that  come 
exclusively  from  the  resonance  we  are  trying  to  select  in  the  same  mass  intervals 
(that  is,  the  true  helicity  angle  distribution  of  the  resonance)  £/0(980),  £/2(i270),  and 
£/4(2050)  respectively.  The  continuum  subtracted  helicity  angle  distribution  of  a 
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resonance  is  being  contaminated  by  the  tails  of  the  other  resonances.  Keeping  this 
in  mind  we  write  2 


C/o(980)  =  */o(980)  +  Ptf2(i270) 

C/2(1270)  -  */o(1270)  +  7*/o(980)  +  <^/4(2050) 


(4.3) 


C/4(2050)  *  fyo(2050)  +  e£/2(i270) 

Where  the  small  numbers  (3,  7,  S,  and  e  are  the  ratios  of  the  number  of  events 
from  a  resonance  in  the  mass  interval  where  the  contamination  is  taking  place  to 
the  number  of  events  from  the  same  resonance  in  the  mass  interval  used  to  select 
it.  Using  the  (mm,  Amm)  values  in  Table  4-2  and  the  fit  in  Figure  4-3  we  obtain 
/3  =  8.5x  1(T3,  7  =  6.0  x  10-2,  5  =  9.7  x  10~2,  and  e  =--  0.13. 

To  obtain  the  background  subtracted  helicity  angular  distributions  we  simply 
invert  the  system  of  equations  expressed  in  Equation  4.3.  The  solution  can  be 
conveniently  expressed  as, 


tf2(1270)  =  r-"^:^(C/o(1270)  -  7C/0(980)  ~  ^(2050) ) 
*/o(980)  =  c/o(980)  -  Ptf2(1270) 


(4.4) 


£/4(2050)  =  C/0(2050)  —  <^/2(1270) 

As  a  check,  if  we  ignore  the  second  order  terms  the  previous  solution  becomes, 


2  Here  the  cross  contamination  between  the  /o(980)  and  the  /4(2050)  is  ignored. 
There  is  no  mathematical  problem  in  including  these  contamination  terms,  but  the 
system  of  equations  would  be  more  complicated  than  it  needs  to  be  since  this  type 
of  cross  contamination  is  negligible. 
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*/a(12TO)  ~  C/0(1270)  -  7C/O(980)  ~  <fc/4(2050) 
<   */o(930)  =  c/o(980)  -  /3C/2(1270) 
*/4(2050)  =  C/o(2050)  -  eC/2(1270) 

which  is  indeed  is  the  solution  when  cross-contamination  is  ignored. 

4.3.3     Statistical  Fit  of  the  Helicity  Angular  Distribution:  J  Assign- 
ments and  their  Probability  Distribution 

For  each  resonance,  we  fit  the  data  to  the  simultaneous  cos  61  and  cos  9n 
helicity  angle  distribution  obtained  from  data  using  Equation  4.4  to  the  helicity 
formalism  prediction,  Equations  23-27.  projected  on  each  angle  and  folded  in 
opposite  directions  around  their  symmetry  axis  in  order  to  show  both  distributions 
on  the  same  plot.  The  corresponding  J  value  of  the  best  fit  for  each  resonance, 
shown  in  Figures  4-8  through  4-10,  is  defined  as  the  J  assignment  (Ja).  We  obtain 
J a  =  1  for  the  /o(980)  which  is  inconsistent  with  the  known  spin  of  the  /o(980) 
which  is  J  =  0.  For  the  /2(1270)  and  the  /4(2050)  we  obtain  Ja  =  2  and  Ja  =  4 
respectively,  which  is  consistent  with  their  known  spins. 

To  have  an  idea  of  how  well  the  angular  distribution  determines  J  among  the 
hypotheses  J  =  0,  1,  2,  3,  4,  we  do  a  statistical  fit  for  each  hypothesis  and  assign 
each  one  a  probability  proportional  to  e-(x2+d.o.f.)/2  where  do  j  are  t^e  degrees 
of  freedom  in  the  fit.  The  resulting  normalized  probability  distributions  give  an 
idea  of  the  assigned  J  significance  and  are  shown  in  Figure  4-11.  In  particular, 
this  figure  shows  that  the  Ja  =  1  for  the  /o(980)  inconsistency  can  be  due  to  a 
statistical  fluctuation. 
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File:  7boh.r/userl/luis/deo3/otp/data_50.ntp 

ID       IDB    Symb              Date/Time  Area 
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Figure  4-2:  Invariant  mass  of  7r+7r-  for  IS  (top)  and  4S  (bottom)  data.  For  the  4S 
data  we  show  a  blow-up  of  the  mass  region  1.5     -     3  GeV/c2  where  the  p*  can  be 


seen. 
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MINUITx2  Fit  to  Plot 

9&3 

UlS->Gamma  h  h.  masspipi  axis 

File:  Generated  internally 
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Plot  Area  Total/Fit             1 544. 1/1542.1 

Fit  status  3 
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Figure  4-3:  Continuum  subtracted  (as  defined  in  Section  3.2.1)  invariant  mass  of 
7T+7r~  from  T(15)  — >  77T+7r~. 
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tile:  Generated  internally 
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Figure  4-4:  Blow  up  of  the  invariant  mass  of  tt+it    in  the  /4(2050)  mass  region, 
including  the  upper  limit  on  the  fj{2220)  shown  as  a  dashed  line. 
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File:/lmp_rnnt/bofif/user1.1uis/c!eo3/>nnJit/u1sgammahh/tertip.tnnp 

ID       IDB    Symb              Date/Time                    Area  Mean  RMS 
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Figure  4-5:  Normalized  probability  distribution  for  different  /j(2220) 
nal  areas.  The  shaded  area  spans  90%  of  the  probability. 
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Figure  4-6:  Contour  plot  of  the  inverse  of  the  average  relative  bin-error  from  back- 
ground subtractions  in  the  /0(980)  (top),  /2(1270)  (middle),  and  /4(2050)  (bottom) 
angular  distribution. 
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MINUITx   Fit  to  Plot  9&82 

U1S->Gamma  h  h  -ABS(C0S(PH0THETA1)  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  606.18  /  606.18 

Func  Area  Total/Fit  436.31  /  436.31 

X2=     25.6  for  20-  1  d.o.f., 

Errors  Parabolic 

Function  1 :  JO 

AREA  218.15  +     42.27     . 

300 


5-AUG-2004  23:06 

Fit  Status  3 

E.O.M.  6.301  E-1 3 
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Figure  4-7:  Angular  distribution  for  the  excess  events  in  the  /o(980)  mass  region. 
The  fit  corresponds  to  J  =  0. 
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MINUIT  %  Fit  to  Plot  9&82 

U1S->Gamma  h  h.  -ABS(COSfPHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  606. 1 3  /  606. 1 8 

Func  Area  Total/Fit  478.89  /  478.89 
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Errors  Parabolic 

Function  1 :  J1 

AREA  239.45  ±     39.31 
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Figure  4-8:  Angular  distribution  for  the  excess  events  in  the  /o(980)  mass  region. 
The  fit  corresponds  to  J  =  1. 
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MINUIT  %   Rt  to  Plot  9&72 

U1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  1673.4 1  1673.4 

Func  Area  Total/Fit  1 733.8  / 1 733.8 

X2=    27.3  for  20-  3d.o.f, 

Errors  Parabolic 

Function  1 :  J2 

AREA  866.97  ±     43.10 
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<I>  0.29838  ±    0.1098 
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Figure  4-9:  Angular  distribution  for  the  excess  events  in  the  /2(1270)  mass  region. 
The  fit  corresponds  to  J  =  2. 
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MINUiTx   Fit  to  Plot  9&92 

U1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  1 59  93  / 1 59.93 

Func  Area  Total/Fit  127.08  / 127.08 

X2=     19.6  for  20-  3d.o.f., 

Errors  Parabolic 

Function  1 :  J4 
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Figure  4-10:  Angular  distribution  for  the  excess  events  in  the  /4(2250)  mass  region. 
The  fit  corresponds  to  J  =  4. 
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Figure  4-11:  Ja  probability  distribution  for  resonances  in  the  /o(980)  (top), 
/2(1270)  (middle),  and  /4(2050)  (bottom)  invariant  mass  region  when  only  the 
hypotheses  Ja  =  0,  1,  2,  3,  4  are  considered. 


CHAPTER  5 
EXCLUSIVE  RADIATIVE  DECAY  T(15)  -*  ^yK+K~ 

In  figure  5-1  the  K+K~  invariant  mass  plot  is  shown  for  both  resonance  and 
continuum  running.  This  figure  also  has  an  inset  showing  the  (JdE/dx{K)  for  both 
tracks  for  events  in  the  1.1  -  3  GeV/c2  region.  This  inset  indicates  that  most 
of  the  events  in  this  mass  region  have  indeed  two  kaons  and  that  the  amount  of 
p  reflection  is  small.  Furthermore,  when  the  ir+n~  invariant  mass  is  plotted  for 
these  events  only  40  out  of  700  can  be  fit  under  a  p  peak  (this  fit  is  what  we  use  to 
calculate  the  pion  faking  kaon  rate  in  Table  3-6). 

Figure  5-2  shows  our  fit  to  the  K+K"  invariant  mass  continuum  subtracted 

plot  as  defined  is  Section  3.2.1  with  the  most  likely  statistical  fit  overlayed  (which 

is  described  in  the  next  section).  The  number  of  events  near  the  (f>  region  (1.01  - 

1.03  GeV/c2)  left  after  the  continuum  subtraction  is  50  ±  70. 

5.1     Statistical  Fit  of  the  Invariant  Mass  Distribution:  Signal  Areas  and 

Their  Significance 

The  results  of  this  section  are  summarized  in  Table  5-1. 

From  the  measurement  of  the  previous  section  we  expect  a  small  contribution 
(«  50  events  assuming  no  interference)  from  /2(1270)  -+  K+K~ .  We  do  find 
some  evidence  for  /2(1270)  -+  K+K~  events  (110  ±  40)  in  the  fit.  We  also  find 
strong  evidence  for  the  resonance  /j(1525)  and  weak  evidence  for  the  /0(1710) 
resonance.  The  /2(1270),  and  /0(1710)  are  fitted  with  their  widths  fixed  to  their 
PDG  values[16]  because  they  have  large  errors  if  allowed  to  float.  The  rest  of 
the  resonances  parameters  are  consistent  with  their  PDG  values[16],  which  are 
ro/5(i525)  =  1525  ±  5  MeV/c2,,  17/(1526)  =  76  ±  10  MeV/c2,  m/o(mo)  -  1715  ± 
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6  MeV/c2,  and  Tfo^n0)  -  125  ±  10  MeV/c2.  We  also  observe  an  excess  of  events  in 
the  2  -  3  GeV/c2  region  which  we  can't  attribute  to  any  known  resonances. 

Significances  of  the  signals  of  the  identified  resonances  in  the  fit  are  calculated 
as  described  in  the  previous  section.  The  T(15)  ™>  7//(2220),  /j(2220)  -»  K+K~ 
upper  limit  is  also  calculated  using  multiple  fits,  except  that  this  time  the  events 
under  the  /j(2220)  are  of  unknown  origin,  so  we  use  a  first  order  polynomial 
allowed  to  float.  The  significance  of  the  excess  of  events  in  the  2  —  3  GeV/c2 
invariant  mass  region  is  calculated  assuming  a  normal  distribution;  we  simply  add 
up  the  number  of  events  in  each  bin  along  with  its  error.  Results  are  shown  in 
Table  5-1. 

Table  5-1:  Results  for  T(15)  -+  jK+R- 


Mode 

Area 

Branching  Fraction  (10  5) 

Significance 

7/2(1270) 

110  ±40 

23  ±8 

5.4  x  10~4(3.3a) 

7/2(1525) 

360+87° 

0  q+0.9 

<  10~45(>  14(j) 

7/o(1710),/o(1710)->^+X- 

75  ±30 

0.35  ±0.14 

7.5  x  10-4(3.2a) 

jK+K~(2-3GeV/c2) 

220  ±  20 

1.03  ±0.12 

8.8a 

7 f j(2220),  f j{2220)  -+  K+K- 

<  10 

<  5  x  10-2 

- 

5.2     Angular  Distribution  of  The  Signal 

In  this  section  we  adapt  the  ideas  presented  in  Section  3.4.3  to  the  T(IS')  — > 
^K+K~  situation. 

The  derived  statistical  errors  from  the  signal  and  background  subtractions  can 
be  used  to  calculate  the  mass  interval  which  best  represents  the  helicity  angular 
distribution.  The  inverse  of  the  expected  average  relative  bin  error  as  a  function  of 
the  mass  interval  is  shown  in  Figure  5-5,  and  the  mass  interval  that  maximizes  it 
are  tabulated  in  Table  5-2. 

The  tails  from  the  resonances  contribute  to  the  continuum  subtracted  helicity 
distributions,  ' 
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Table  5-2:  Mean  masses,  widths  inmass  and  inverse  of  the  average  bin  relative 
error  (F)  from  background  subtractions  for  the  angular  distribution  of  different 
resonances.  Standard  mean  masses  and  widths,  corresponding  to  1I\  are  taken 
from  the  fit  in  Figure  5-2,  and  are  labeled  with  the  subscript  "s" ,  while  those  that 
maximize  F  are  labeled  with  the  subscript  "m" .  The  last  column  shows  the  factor 
by  which  the  effective  statistics  increase. 


Resonance 

.... 
(ms,  Am,) 

F(ras,Aras) 

(ram,  A??im) 

F(mm,  Amm) 

F2(ms,Ama) 

/2(1270) 
£(1525) 

/o(1710) 

(1.276,0.185) 
(1.540,0.085) 
(1.760,0,125) 

3.5 
9.6 
3.1 

(1.300,0.100) 
(1.565,0.100) 
(1.780,0.095) 

4.1 
9.8 
3.3 

1.37 
1.04 
1.13 

(5.1) 


C/2(1270)  -  */a(1270)  +  /^/2(1525) 

C/2(1525)  ~  */£(1525)  +  7*/2(1270)  +  ^*/o(1710) 

cjb(mo)  =  tfo{mo)  +  ei/2(i525) 

Where  again  the  small  numbers  /?,  7,  5,  and  e  are  the  ratios  of  the  number 
of  events  from  a  resonance  in  the  mass  interval  where  the  contamination  is  taking 
place  to  the  number  of  events  from  the  same  resonance  in  the  mass  interval  used 
to  select  it.  Using  the  (mm,  Amm)  values  in  Table  5-2  and  the  fit  in  Figure  5-2  we 
obtain  /?  =  1.3  x  10"2,  7  -  0.14,  5  =  0.11,  and  e  =  0.19. 

The  background  subtracted  helicity  angular  distributions  are, 


(5.2) 


%1525)  =   1-/^-^(^(1525)  -  7C/2(1270)  ~  <5c/o(m0)) 
i/2(1270)  =  C/2(i270)  -  Pt/^1525) 
</o(1710)  =  c/o(1710)  -  ^/2(1525) 

The  the  best  fit  for  each  resonance  and  the  excess  of  events  in  the  2—3  GeV/c2 
are  shown  in  Figures  5-6  through  5-9.  The  best  spin  assignment  for  the  /2(1270)  is 
Ja  =  2,  for  the  /2(1525)  it  is  Ja  =  2,  for  the  it  is  /0(1710),  and  it  is  Ja  =  1  for  the 
excess  of  events  in  the  2  —  3  GeV/c2  mass  region.  The  Ja  =  2  value  for  the  /o(1710) 
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is  inconsistent  with  its  known  spin.  Also,  examination  of  the  normalized  helicity 
amplitudes  for  the  /2(1270)  reveals  that  they  are  inconsistent  with  those  obtained 
for  the  /2(1270)  in  the  7r+7r~  mode. 

The  assigned  J  probability  distributions  are  shown  in  Figure  5-11.  They  reveal 
that  the  inconsistencies  in  the  /2(1270)  and  /0(1710)  are  not  significant  and  can  be 
attributed  to  the  statistical  uncertainty. 


Fiie:7bohr/user1/luis/deo3/ntp/data_50.ntp 

ID       IOB    Symb             Date/Time 

Area 

Mean 

R.MS. 

1           1           1          040807/0219 

1274. 

1.320 

04283 

1            2            1           040807/0219 

2031. 

1.156 

0.2585 

1          10            1           040807/0235 

705.0 

0.2246 

1.208 

0.1480 

1.170 

! 

© 


1 
p 
> 


600 


400 


200 


1200 


800 


400 


-t — i — i — r 


t — i — r—r 


0.9 


|     II    1     I  -|"l  -1  ■  1  -1     | 


6 

h 

-6 


■on 


1.4 


1.9 
Mass  (GeV/c2) 


2.4 


2.9 


Figure  5-1:  Invariant  mass  of  K+K~  for  IS  (top)  and  4S  (bottom)  data.  For 
the  4S  data  the  inset  shows  the  <JdE/dxi_K)  fc>r  both  tracks  for  events  in  the 
1.1    —    3  GeV/c2  mass  region.  This  inset  indicates  that  most  of  the  events  are  con- 
stant with  having  two  kaons. 
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Figure  5-2:  Continuum  subtracted  invariant  mass  of  K+K~ 
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Figure  5-3:  Blow  up  of  the  /j(2220)  region.,  with  the 
laid.  The  mass  and  width  are  taken  to  be  ra/^2220) 
r/j(2220)  =  17  MeV/c2  as  in  [17]. 
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File:  /tmp_ mnt/bohr/userl/luis/cleo3/mn  fitAilsgammahh/temp.tmp 

ID       IDB    Symb              Date/Time                    Area  Mean  H.M.S. 

1          29            1           040807/1426                    1.000  4.746  3.582 

1          ?5            1           040807/1426                  0.9054  3.971  2.723 
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Figure  5-4:  Normalized  probability  distribution  for  different  /j(2220) 
signal  areas.  The  shaded  area  spans  90%  of  the  probability. 
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Figure  5-5:  Contour  plot  of  the  inverse  of  the  average  relative  bin-error  from  back- 
ground subtractions  in  the  /2(1270)  (top),  /2(1525)  (middle),  and  /0(1710)  (bot- 
tom) angular  distributions. 
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MINUITx   FittoPlot  9&82 

U1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  517.49/517.49 

Func  Area  Total/Fit  129  01  / 129.01 


%=     17.8  for  20-  3d.o.f., 


5-AUG-2004  23:16 

Fit  Status  3 

E.D.M.  5.475F.-Q5 

C.L=  40.4% 


Errors 
Function  1 :  J2 
AREA 

e 
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100 


64.505 
7.95385E-03 
1.4108 
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Figure  5-6:  Background  subtracted  K+K    angular-  distribution  in  the  /2(1270) 
mass  region- as  denned  in  the  text.  The  fit  corresponds  to  J  =  2. 
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MINUIT  x   Fit  to  Plot  9&72 

U1S->Gamma  h  h,  -ABS(COS(PHOTHETA))  axis 
Pile:  Genefaled  internally 
Plot  Area  Total/Fit  284  44  /  284.44 

Func  Area  Total/Fit  367.87  /  367.87 

X2=    20.9  for  20-  3d.o.f., 


5-AUG-20O4  23:16 

Fit  Status  3 

E.D.M.  3.771  £-06 

C.L.=  23.3% 


Errors 
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Figure  5-7:  Background  subtracted  K+  K    angular  distribution  in  the  f2{\525) 
mass  region  as  denned  in  the  text.  The  fit  corresponds  to  J  =  2. 
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MINUIT  y.   Fit  to  Plot  9&92 

iJ1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  1 29.47  / 1 29.47 

Furic  Area  Total/Fit  79.031  /  79.031 

X2=    28.0  for  20-  3d.o.f., 

Errors  Parabolic 

Function  1 :  J2 

AREA  39.516  ±     10.80 

e  2.99860E-04  ±    0.3847 

<S>  1.5696  ±    0.2633 

50 
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Fit  Status  3 
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Figure  5-8:  Background  subtracted  K+K    angular  distribution  in  the  /0(1710) 
mass  region  as  defined  in  the  text.  The  fit  corresponds  to  J  =  2. 
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MINUITx z  Fit  to  Plot  9&102 

U1S->Gammah  h.  ABS(GOS(PHOTHErrA))  axis 
f-'ile:  Generated  internally 
Plot  Area  Total/Rt  41 1 .73  /  41 1 .73 

Fjnc  Area  Total/Fit  370.86  /  370.86 

X2=     18.9  for  20-  2d.o.f., 

Errors  Parabolic 

Function  1:J1 

AREA  165.43  ±     17.19 

0  0155<t2  ±   0.3002 
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Figure  5-9:  K+K~  angular  distribution  for  events  within  the  2  -  3  GeV/c2 
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MINUIT  x   Fit  to  Plot  9&92 

U1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  1 29.47  / 1 29.47 

Func  Area  Total/Fit  62. 1 94  /  62. 1 94 

%2=     33.3  for  20-  1  d.o.f., 

Errors  Parabolic 

Function  1 :  JO 

AREA  31.097  +     10.92 
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Figure  5-10:  Background  subtracted  K+K~  angular  distribution  in  the  /0(1710) 
mass  region  as  defined  in  the  text.  The  fit  corresponds  to  J  =  0. 
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Figure  5-11:  Ja  probability  distribution  for  resonances  in  the  fi{  1270) (top), 
/2(1525)  (middle-top),  /o(1710)  (middle-bottom)  invariant  mass  region  and  the 
excess  of  events  in  the  2  —  3  GeV/c2  (bottom)  region  when  only  the  hypotheses  Ja 
—  0,  1,  2,  3,  4  are  considered. 


CHAPTER  6 
EXCLUSIVE  RADIATIVE  DECAY  T(15)  -» 1Pp 

In  figure  6-1  the  pp  invariant  mass  plot  is  shown  for  both  IS  and  4S  data,  with 
an  inset  showing  that  most  of  the  events  indeed  have  a  proton  and  an  anti-proton 
and  that  the  p  and  4>  reflections  are  small.  The  enhancement  at  3.1   GeV/c2  in 
the  4S  plot  corresponds  to  the  process  e+e"  — >  ^J/tp  with  J/ip  -*  pp.  This 
enhancement  is  not  as  pronounced  in  the  IS  plot  because  the  IS  data  only  has 
32.4%  of  the  luminosity  of  the  4S  data.  The  number  of  events  in  the  J/tp  invariant 
mass  region  after  continuum  subtraction  is  6  ±  3  and  is  consistent  with  0. 

Figure  6-2  shows  the  continuum  subtracted  invariant  pp  mass  distribution 
(as  defined  in  Section  3.2.1)  with  a  90%  confidence  level  upper  limit  for  /j(2220) 
overlaid.  A  direct  fit  to  the  fj  (2220)  yields  12  ±  5  events. 

There  is  an  excess  of  events  in  the  continuum  subtracted  invariant  mass  plot. 
We  measure  this  excess  and  the  upper  limit  of  e+e~  -»  7/j(2220),  /j(2220)  ->  pp 
the  same  way  we  measured  the  excess  of  events  inside  2  GeV/c2  <  rn(K+K~)  < 
3  GeV/c2  region  and  the  /j(2220)  upper  limit  in  Section  5.1.  Results  are  shown  in 
Table  6-1. 

The  pp  angular  distribution  for  the  mass  range  2  GeV/c2<  m(pp)  <  3  GeV/c2 
is  shown  in  Figure  6-4  with  the  most  likely  J  assignment,  Ja  =  1,  fit  overlayed. 
The  probability  distribution  for  Ja  is  shown  in  Figure  6-5. 

Table  6-1:  Results  for  T(15)  -*  7pp. 


Mode 

Area 

B.F.  or  90  %  U.L.  (lO"5) 

Significance 

7PP  (2-3  GeV/c2) 
7/j(2220),/j(2220)-+pp 

85  ±18 
<20 

0.41  ±0.08 
<  9  x  10-2 
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Rle:7bohr/user1/luis/deo3/ntp/data_50.ntp 
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Figure  6-1:  Invariant  mass  of  pp  for  IS  (top)  and  4S  (bottom)  running.  For  the  4S 
running  the  inset  is  consistent  with  the  events  having  a  proton  and  an  anti-proton. 
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File:  Generated  internally 
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Figure  6-2:  Invariant  mass  of  pp.  The  plot  is  continuum  subtracted  and  efficiency 
corrected.  An  overlay  with  the  90%  confidence  level  upper  limit  for  fj (2220) 
is  shown.  The  mass  and  width  are  taken  to  be  777.^(2220)        =        2.234  GeV/c2  and 
r/j(2220)  =  17  MeV/c2  as  in  [17]. 
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File:/tmp_mnt/bohr/user1/luis/deo3/mn_fit/u1sgammahh/temp.tmp 

ID       IDB    Symb              Date/Time                    Area  Mean  R.M.S. 

1          29            1           040807/1446                    1.000  12.78  5.430 

1           25             1            040807/1446                    0.9110  11.80  4.590 


0.060 


§ 

X> 

'B   0.040 

1/3 

-3 


■8 

2 

cu 


0.020 


0.000 


t 1 1 r 


T 1 ! P 


25 


50 


Area 


Figure  6-3:  Normalized  probability  distribution  for  different  /jr(2220)    — ♦    pp  signal 
areas.  The  shaded  area  spans  90%  of  the  probability. 
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MINUITjc   FittoPlot  9&22 

U1S->Gamma  h  h.  -ABS(COS(PHOTHETA))  axis 
File:  Generated  internally 
Plot  Area  Total/Fit  256.69/256.69 

Func  Area  Total/Fit  1 34.57  / 1 34.57 

%2=     19.5  for  20-  2d.o.f., 

Errors  Parabolic 

Function  1:J1 

AREA  67.287  +     10.91 

0  1.20611E-03  +    0.4148 
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Figure  6-4:  pp  angular  distribution  for  the  excess  of  events  in  the  mass  range 
2  GeV/c2<  m(pp)  <  3  GeV/c2. 
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Figure  6-5:  Ja  probability  distribution  for  the  excess  of  events  in  the  mass  range 
2  GeV/c2<     m(pp)     <     3  GeV/c2  when  only  the  hypotheses  Ja  =  0,  1.  2,  3,  4  are 
considered. 


CHAPTER  7 
SYSTEMATIC  UNCERTAINTIES 

Systematic  uncertainties  are  any  sources  of  experimental  uncertainty  other 
than  the  statistical  ones.  Limits  on  the  accuracy  of  our  detector  simulation  and  any 
physical  processes  that  interferes  with  the  experimental  measurement  are  typical 
examples  of  systematic  uncertainties.  These  uncertainties  need  to  be  identified, 
quantified,  and  if  possible,  corrected. 

7.1     Cuts 

From  now  on  we  report  individual  cut  efficiencies  relative  to  the  events  that 
survive  all  the  other  cuts. 

According  to  our  MC,  most  of  the  skim  cuts,  except  for  the  hardGam  re- 
quirement, are  nearly  100%  efficient.  Therefore,  such  cuts  should  not  be  a  source 
of  systematic  uncertainty.  The  only  skim  cut  worth  taking  a  closer  look  at  is  the 
hardGam  cut,  which  is  about  93%  efficient  in  MC.  Measuring  this  efficiency  in  our 
data  is  not  possible  because  events  not  classified  as  hardGam  are  not  in  the  data 
to  begin  with.  Closer  examination  of  our  MC  reveals  that  that  most  of  the  7% 
inefficiency  in  hardGam  comes  from  the  eOverPl  cut  (4%),  some  from  the  Sh2  cut 
(2%),  and  the  rest  (1%)  from  the  other  cuts  present  in  hardGam.  We  can  measure 
the  efficiency  of  eOverPl  in  data  by  looking  at  the  eOverP2  distribution  of  data 
tracks  from  p  and  <p  decay.  We  are  satisfied  by  this  check  on  the  MC  modeling  (see 
Tables   7-1  and   7-2),  and  we  won't  measure  how  well  MC  models  the  rest  of  the 
hardGam  cuts,  which  are  97%  efficient  in  MC. 
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To  quantify  the  quality  in  the  MC  modeling  of  our  analysis  cuts  in  the  ix 
and  K  modes  we  use  the  p  and  <j)  signals  present  in  our  data1  .  In  this  case  the 
process  is  straightforward.  We  measure  the  p  and  <p  signal  signals  in  data  and  MC 
over  a  floating  background  function  with  all  cuts  in  place,  and  with  all  in  place 
cuts  except  the  one  under  consideration.  Prom  these  numbers  we  calculate  the 
effective  efficiency  of  the  cut.  The  differences  between  data  and  MC  are  taken  as 
the  systematic  errors  which  are  added  in  quadrature.  Results  are  shown  in  Tables. 
7-1-  7-2.  The  4-momentum  cut  does  not  appear  because  its  efficiency  is  close  to 
100%. 


Cuts 


MCeff. 


Data  eff. 


Systematic  Error 


eOverP2 

QED  e+e"  — ►  7e+e~  suppresion 

QED  p+p~  — >  7P+p~  suppresion 

Hadron  separation 


96.2 
93.1 
75.3 

96.9 


96.6 
94.2 

77.7 
96.8 


0.4 
1.2 
3.2 
-0.1 


Overall  analysis  cut  systematic  error 


±3 


Table  7-1:  Efficiencies  for  T(15)      — ►      ~/n+n~  for  flat  signal  MC,  efficiencies  from 
data  (p),  and  the  derived  systematic  error  in  %.  Efficiencies  are  reported  as  the 
number  of  signal  events  after  all  cuts  divided  by  the  number  of  signal  events  with 
all  cuts  except  the  one  under  consideration.  Statistical!  errors  are  0.1%  or  less. 


For  the  proton  case  we  don't  have  a  clean  sample  with  high  statistics  of 

7PP  events  in  data.  By  extension  we  take  the  systematic  error  in  this  mode 


e+e 

to  be  10%. 

7.1.1     Justification  of  the  DPTHMU  Cut 

The  reason  we  prefer  to  use  DPTHMU  <  5,  instead  of  the  more  traditional 
(see   [18])  DPTHMU  <  3  used  in  CLEO  II,  is  that,  for  some  unknown  reason,  our 
CLEO  III  MC  has  too  many  -k  tracks  with  3  <  DPTHMU  <  5. 


1  A  study  using  Kg  or  A  signals  from  hadronic  environments  would  have  larger 
statistics,  but  is  problematic  because  the  large  number  of  tracks  and  showers  artifi- 
cially decrease  the  cut  efficiency.  See  [28]  Appendix  A. 3  for  an  example. 
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Cuts 

MCEff. 

Data  Eff. 

Systematic  Error 

eOverP2 

QED  e+e~  — »  je+e~  suppresion 

QED  p+p~  — ►  ,yp+p~  suppresion 

Hadron  separation 

98.5 
98.1 
93.5 
88.4 

100 
99.3 
96.4 
82.2 

1.5 
1.3 

3.1 

-7.0 

Overall  analysis  cut  systematic  error 

±8 

Table  7-2:  Efficiencies  for  T(15)     -»     ^K+K~  for  flat  signal  MC,  efficiencies  from 
data  ((/)),  and  the  derived  systematic  error  in  %.  Efficiencies  are  reported  as  the 
number  of  signal  events  after  all  cuts  divided  by  the  number  of  signal  events  with 
all  cuts  except  the  one  under  consideration.  Statistical!  errors  are  0.1%  or  less. 


To  observe  this  fact  we  first  select  a  relatively  clean  sample  of  pion  tracks 
by  requiring  the  event  to  have  a  7r+7r~  invariant  mass  consistent  with  the  p  mass, 
and  to  pass  all  our  analysis  cuts  except  for  the  cut  on  DPTHMU  on  one  track. 
For  such  events  we  plot  the  rate  as  a  function  of  momentum,  at  which  the  n  track 
whose  DPTHMU  cut  we  released  has  DPTHMU  >  3  and  DPTHMU  >  5.  To 
increase  our  statistics  we  do  this  procedure  twice,  once  for  each  track,  and  average 
the  fake  rate.  Figure  7-1  shows  the  results.  Clearly  the  CLEO  III  MC  we  are  using 
has  some  problem  modeling  the  DPTHMU  <  3  cut. 

To  keep  the  systematic  error  low  we  choose  a  cut  at  DPTHMU  <  5 
(DPTHMU  <  3  gives  a  systematic  error  of  about  20%).  This  does  not  change  the 
efficiency  in  data  very  much,  but  it  increases  the  efficiency  i€^3?ted  by  MC,  bring- 
ing it  closer  to  reality.  The  increase  in  p  fakes  after  loosening  the  cut  is  estimated 
to  be  low  using  QED  MC  (see  Figure  7-2). 

7.2     Angular  Distribution  of  Signal 

The  photon  and  tracks  from  the  process  T(IS')  — »  7X  with  X  -*  h+h~  have 
a  different  angular  distribution  than  that  of  flat  MC.  Examples  of  possible  angular 
distributions  are  shown  in  the  appendix. 
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Figure  7-1:  Pion  faking  muon  fake  rates  for  a  cut  on  DPTHMU      <     3  (top)  and 
DPT  HMD  <  5  (bottom).  MC  is  shown  as  solid  circles  while  p  from  data  is  shown 
as  hollow  circles.  Fake  rates  are  reported  relative  to  all  events  that  pass  all  cuts, 
except  for  the  DPTHMU  cut  for  one  of  the  tracks. 

Figure   7-3  shows  the  efficiency  in  flat  MC  as  a  function  of  cos  #7  and 
cos9h+2  ■  Note  that  the  K  and  p  modes  are  nearly  insensitive  to  the  track  an- 
gular distribution,  while  the  it  mode  is  more  sensitive  to  6^+  ■  This  happens  because 
of  the  stronger  muon  rejection  cut  in  the  7r  mode. 

We  measure  the  systematic  effects  of  flat  MC  efficiency  by  convoluting  each 
plot  in  Figure  7-3  with  different  possible  angular  distributions  calculated  in  the 


2  61  and  $h+  are  the  helicity  angles  of  the  sequential  decay,  defined  in  the  ap- 
pendix (see  Figure  2). 
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Figure  7-2:  Muon  faking  pion  fake  rate  (in  %)  for  DPTHMU    <    3  (solid  squares) 
and  DPTHMU  <  5  (solid  circles).  There  is  a  complete  overlap  in  the  last  bin. 

appendix.  Tables  7-3-7-4  show  the  necessary  correction  factors  relative  to  flat  MC 
efficiency  for  decays  with  definite  7  and  X  helicities  A7  and  Ax,  due  to  the  non-flat 
photon  and  hadron  distributions.  We  call  these  factors  cxxx   (j)  and  exxXx(h) 
respectively. 

The  efficiency  of  a  decay  with  definite  A7,  Xx  can  be  obtained  using  the  flat 
MC  efficiency  corrected  by  a  factor  e\Xx  —  £\x\x(,l)  x  exxxx{h)-  In  general,  the 
final  state  is  a  mixture  of  all  possible  A7;  A^  pairs,  and  the  efficiency  correction 
factor  is, 


eJx  =  cos2  0  cos2  $eft  +  sin2  0efr  +  cos2  0  sin2  toe!? . 


-10 


11 


(7.1) 
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Figure  7-3:  Flat  MC  efficiency  as  a  function  of  cosO^  (left  column)  and  as  a  func- 
tion of  cos6h+  (right  column).  Top  row  corresponds  to  the  pion  mode,  middle  row 
to  the  kaon  mode  and  bottom  row  to  the  proton  mode. 

The  fits  in  Figures  4-8  and  4-10  and  5-6-5-8  measure  the  pair  (0,  $).  These 
values  are  summarized  in  Table  7-5,  where  they  are  used  to  obtain  eJx  for  each 
mode. 

The  pair  (0,  $)  carry  an  error  which  is  a  source  of  systematic  uncertainty. 
We  calculate  this  systematic  uncertainty  by  inspecting  the  differences  in  efficiency 
when  (0, 4>)  move  away  from  the  value  which  gives  the  minimum  chi- squared,  Xmini 
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Mode 

aw  or  an 
9  =  0 

an 

9  =  - 

W          2 

K 
P 

0.910 
0.925 
0.922 

1.18 
1.15 
1.16 

Table  7-3:  e^A   (7),  efficienciy  correction  with  respect  to  flat  MC  factors  due  to 
the  non-flat  angular  distribution  of  the  photon  in  T(1S)  -»  7X  for  different  X  spin 
values  (Jx)  and  7,  X  helicities  (A7,Ax)- 


Mode 

■/* 

«10 

9  =  0,  *  =  0 

an 

9  =  - 
w      2 

812 

9  =  0,  $  =J 

0 

1.00 

- 

- 

1 

0.785 

1.11 

- 

X  -+  tt+tt- 

2 

0.829 

0.934 

1.15 

3 

0.898 

0.853 

1.03 

4 

0.899 

0.863 

0.922 

0 

1.00 

- 

- 

x  -» #+#- 

1 
2 

0.944 
0.950 

1.03 
0.988 

1.04 

3 

0.957 

0.971 

1.01 

4 

0.965 

0.962 

0.992 

Table  7-4:  exxXx(h),  efficienciy  correction  factors  with  respect  to  non-flat  MC  due 
to  the  track  angular  distribution  in  T(15)  — >  7X,  with  X  — »  h+h~  for  different  X 
spin  values  (Jx)  and  7,  X  helicities  (A7,  Ax)- 


under  the  condition  x2  <  Xmm  +  1  3  •  F°r  Jx  >  1  both  (9,  $)  are  free  to  move, 
defining  a  surface  in  the  (9,  $)  plane.  These  surfaces  are  shown  in  Figure  7-4  for 
the  /2(1270),  /4(2050),  and  £(1525). 

At  this  point  we  can  check  whether  (9,  $)  depend  on  the  mass  of  the  de- 
cay. We  split  the  /2(1270)  and  the  /2(1525)  into  a  high  mass  and  a  low  mass 
region.  The  plots  of  the  error  surfaces  of  the  measured  (9,  $)  show  no  significant 
separation  for  the  different  mass  regions  (see  Figures  7-5  and  7-6). 


3  In  a  two  dimensional  linear  problem  such  a  set  of  points  defines  the  surface  of 
the  standard  error  ellipse. 
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Results  for  the  correction  factor  and  its  systematic  error  are  shown  in  Table  7- 
5. 

Upper  limits  on  T(15)  -»  7/2(2220),  /2(2220)  -»  h+h~  are  changed  to 
include  the  angular  distribution's  effect  on  efficiency.  Since  we  can't  measure 
the  helicity  amplitudes  in  this  case,  we  choose  the  worst  possible  case  where  the 
corrected  efficiency  is  lowest.  This  always  corresponds  to  0  =  $  =  0.  Results  are  in 
Table  7-6. 
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Figure  7-4:  Surfaces  in  the  (0,  $)  plane  used  to  determine  the  systematic  uncer- 
tainty in  the  efficiency  correction  factor  for  the  .modes  with  /a  (1270),  /4(2050), 
and  /2(1525).  The  /4(2050)  surface  may  seem  large,  but  when  drawn  in  spherical 
coordinates  it  is  a  small  "north  pole  cap" . 
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Figure  7-5:  Measured  (6,$)  surfaces  for  the  /2(1270)  high  and  low  mass  regions. 
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Figure  7-6:  Measuret.  (0,$)  surfaces  for  the  /2(1525)  high  and  low  mass  regions. 
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7.3     Different  Hadronic  Fake  Rates  Between  IS  and  4S 

Differences  in  the  fake  rates  of  hadronic  events  between  the  IS  and  45 
data  are  another  potential  source  of  systematic  errors  when  doing  a  continuum 
subtraction. 

In  principle  this  is  a  second  order  effect  since  the  contamination  is  proportional 
to  the  difference  of  two  small  numbers.  However  since  we  have  such  a  large  number 
of  p  and  (j)  events  in  the  underlying  IS  continuum,  we  should  quantify  any  possible 
signal  contamination.  Section  3.2.2  outlines  the  three  cases  we  need  to  worry 
about;  p  events  contaminating  the  K+K~  invariant  mass  plot,  p  contaminating  the 
pp  invariant  mass  plot,  and  cp  events  contaminating  the  pp  invariant  mass  plot. 

It  is  reasonable  to  assume  that  the  IS  and  4S  fake  rates  stay  within  50% 
of  each  other  (see  Section  3.3.2  for  some  examples  from  MC  measurements). 
With  this  assumption,  taking  np  «  20000  and  n^  «  2000  in  the  IS  data,  and 
the  fake  rates  measured  in  data  (see  Table  3-6)  we  can  calculate  the  systematic 
uncertainties.  We  expect  ±27  events  from  p  to  contaminate  the  K+K~  invariant 
mass  plot.  However,  these  events  can  be  ignored  as  source  of  systemtic  uncertnty 
because  90%  of  them  fall  in  the  1  -  1.5  GeV/c2  region,  below  the  /2(1525)  peak. 
We  also  expect  ±7  from  p  and  ±3  events  from  4>  contaminating  the  pp  invariant 
mass  plot.  This  represents  a  ±8%  systematic  error  for  the  T(15)  — *  'ypp  mode. 

7.4     Other  Systematic  Sources 

Besides  the  systematic  uncertainties  described  above,  we  also  add  a  2% 
systematic  effect  from  track  finding  (1%  per  track)  and  5%  from  the  number  of 
T(15). 

We  find  no  evidence  for  non-resonant  hadron  pairs  that  could  interfere  with 
the  resonances.  There  is  a  possible  interference  between  /2(1275)  and  ^(1525) 
because  they  both  are  in  the  J  =  2  state.  From  the  /2(1275)  measurement  in  the 
7r+7r~  mode  we  expect  3  events  from  this  source  in  the  /2(1525)  mass  region  1.45  — 


98 

1.6   GeV/c2.  In  the  most  extreme  cases  we  have,  N  <x  1 1^/(1525)  1 12  +  I IAfe(i270) 1 12  ± 
2||-^/'(i525)llll-^/2(i270)||j  where  A.'(1525)  and  Af^mo)  are  the  amplitudes  associated 
with  /2(1525)  and  /2(1270)  respectively.  From  the  data  we  have  iV  ss  300,  and 
using  the  fact  that  ||v4/2(i270)!|2  contributes  with  3  events,  the  contribution  from 
11^/'      5)H2  can  range  from  400  to  200  events.  This  represents  a  30%  variation  from 
the  central  value.  Since  this  is  the  extreme  interference  case,  we  take  half  of  this 
value,  15%,  as  a  reasonable  systematic  uncertainty  from  possible  interference  in  the 
£(1525). 

The  possible  decay  Y(IS')  — >  pir  could  contaminate  our  n+  ir~  invariant  mass 
plot,  since  a  n°  with  a  momentum  greater  than  4  GeV/c2  looks  very  much  like  a 
photon.  Using  phase  space  MC  and  isospin  symmetry,  we  estimate  that  about  40% 
of  possible  T(15)  — ♦  pn  have  such  a  n°.  The  latest  measured  upper  limit    [24] 
for  this  decay  mode  is  1.9  x  10"~6.  This  gives  an  upper  limit  of  16  events  in  the 
efficiency  corrected  ir+  tC  invariant  mass  plot,  of  which  10  events  fall  in  the  p  mass 
region,  behind  the  /o(980)  mass  region.  The  remaining  potential  6  events  could 
contaminate  the  rest  of  the  invariant  mass  region.  None  of  the  signal  regions  is 
affected  by  more  than  1%. 

To  test  the  robustness  of  the  statistical  fits,  we  redo  them  with  all  masses  and 
widths  of  possible  resonances  required  to  stay  within  1<j  of  their  nominal  PDG 
values   [16].  The  results  from  fitting  the  data  in  this  way  are  within  the  statistical 
errors  when  compared  to  the  fit  results  summarized  in  Table  4-1  and   5-1.  We  do 
not  observe  a  systematic  difference  between  these  two  fitting  techniques. 

The  branching  ratios  for  /2(1270)  -*  tt+tt",  /4(2050)  ->  tt+tt-,  and  /a (1525)  -> 
K+K~  are  taken  from  the  PDG  [16]  and  contribute  to  the  systematic  uncertainty 
with  +3,  ±9,  and  ±3  percent  respectively. 


99 


7.5     Overall  Systematic  Uncertainties 

We  combine  the  systematic  errors  of  each  section  in  quadrature.  Tables  7-7- 

7-8  summarize  the  systematic  errors. 

Table  7-7:  Systematic  errors  expressed  as  %  from  the  different  sources  described 
in  the  text.  Not  shown  in  the  table,  but  included  in  the  total,  are  a  systematic 
uncertainty  of  5%  from  the  number  of  Y(1S),  a  2%>  systematic  uncertanty  from 
MC  tracking  in  all  the  modes,  a  15%  systematic  uncertainty  in  the  7/^(1525)  from 
possible  interference  with  7/2  (1270),  an  8%  systematic  uncertainty  in  the  7pp 
mode  from  possible  hadronic  contamination,  and  the  systematic  uncertanties  in  the 
/2(1270),  /4(2050),  and  /2(1525)  hadronic  branching  fractions. 


Mode 

Analysis 
Cuts 

Angular 
Distribution 

Total 

7/o(980) 

7/2(1270) 

7/4(2050) 

±3 
±3 
±3 

±3 

+2 

±6 

±7 
±7 

7/2(1525) 

7/o(1710) 

7.K+A'-(2-3GeV/c2) 

±8 
±8 
±8 

±1 

0 

+1 
-3 

±20 
±10 
±10 

7Pp(2  -  3  GeV/c2) 

±10 

-3 

±14 

Table  7-8:  Increase  in  the  upper  limits  of  T(15)  -»  7/j(2220)  for  different  fj(2220) 
decay  modes  due  to  the  possible  systematic  effects  described  in  the  text  added  in 
quadrature.  Not  shown  in  the  table,  but  included  in  the  total,  are  a  5%  contribu- 
tion from  the  number  of  T(IS')  and  of  2%  contribution  from  MC  tracking. 


Mode 

Analysis 
Cuts 

Angular 
Distribution 

Total 

fj(2220)  -»  tt+tt- 

/j(2220)  -»  K+K~ 

f  j  (2220)  ->  pp 

±3% 
±8% 
±10% 

+33% 
+14% 
+9% 

+34% 
+17% 

+14% 

CHAPTER  8 
RESULTS  AND  CONCLUSION 

We  report  on  a  new  search  for  two-body  radiative  T(15)  decays.  We  place 
stringent  limits  on  the  production  of  the  /j(2220)  particle  in  the  pion  and  kaon 
modes,  and  a  less  stringent  limit  in  the  proton  mode  where  some  excess  of  events  is 
observed  in  the  region  of  interest. 

In  the  decay  channel  T(IS')  — >  77r+7r~  we  find  clear  evidence  for  the  resonance 
/2(1270)  and  measure  a  branching  fraction  of  (13.3  ±  1.4)  x  1CT5,  which  is 
consistent  with  the  earlier  CLEO  measurement  of  (7AtH)  x  10-5  [18].  The  angular 
distributions  of  the  photon  and  tracks  strongly  indicate  that  the  hadron  pairs  we 
assign  to  the  /2(1270)  are  indeed  in  a  J  =  2  state,  and  that  nature  prefers  to 
produce  the  /2(1270)  with  0  helicity.  In  contrast,  for  J/ip  ->  7/2 (1270)  it  was 
found  [29]  that  the  /2(1270)  is  produced  at  equal  rates  with  both  helicity  0  and 
helicity  1,  and  at  a  rate  consistent  with  0  for  helicity  2,  but  more  recently  [30] 
measured  helicity  0  dominance  for  this  same  mode. 

There  is  a  barely  significant  (4.3a)  excess  of  events  in  the  /o(980)  — ►  7r+7r~ 
invariant  mass  region.  The  angular  distributions  of  the  photon  and  tracks  indicate 
that  this  excess  of  events  is  in  a  J  =  1  state,  rather  than  J  —  0  which  would  be  the 
case  if  the  excess  were  due  to  /o(980)  decays.  It  may  be  that  this  excess  is  due  to 
non-resonant  77r+7r~  production.  We  conclude  that  more  data  is  needed  to  resolve 
this  situation. 

We  find  weak  evidence  for  the  production  of  the  resonance  ^(2050)  in  the 
T(l£)  — »  77r+7r-  decay  channel. 

In  the  decay  channel  T(IS')  — >  ryK+K~  we  find  strong  evidence  for  the 
production  of  the  resonance  /2(1525)  and  weak  evidence  for  the  production  of 
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/o(1710).  The  photon  and  track  angular  distributions  show  that  the  two  tracks 
attributed  to  the  /2(1525)  are  indeed  in  a  3  «=  2  state  and  that  the  /2(1525)  is 
produced  mostly  with  helicity  0. 

There  is  also  evidence  of  an  excess  of  events  in  the  2  —  3  GeV/c2  K+K~ 
invariant  mass  region,  which  we  cannot  attribute  to  any  known  resonances. 

Finally,  we  also  find  evidence  of  an  excess  of  events  in  the  channel  T(15)  — > 
7PP,  which  we  cannot  attribute  to  any  known  resonances.  We  find  no  evidence 
of  an  enhancement  near  pp  threshold.  Such  an  enhancement  was  recently  re- 
ported by  [31]  in  the  3  jib  system,  and  is  currently  being  interpreted  (See 
e.g.  [32],  [33],  [34]). 

When  comparing  the  radiative  T(15)  decays  studied  in  this  analysis  to  ra- 
diative 3 1 ib  decays  we  observe  a  suppression  ratio  of  0.09  that  is  in  reasonable 
agreement  with  with  the  naive  expectation  of  0.04  obtained  from  scaling  argu- 
ments. 

Table  8-1  summarizes  our  results. 

While  we  find  hadronic  resonances  at  a  level  consistent  with  the  estimate 
of  Section  1.2,  we  do  not  observe  any  glueball  states.  This  is  not  the  result  we 
would  expect  from  a  naive  interpretation  of  QCD.  as  the  two  gluons  from  a  T(1S) 
radiative  decay  can  form  a  bound  state  directly  and  need  to  go  through  at  least 
two  strong  interaction  vertexes  to  form  a  hadron. 

Internal  consistency  in  QCD  predicts  the  existence  of  glueballs.  Even  though 
theoretical  calculations  of  the  glueball  mass  spectrum  exist,  a  clear  glueball 
observation  has  not  been  made  yet.  We  conclude  that  more  glueball  searches,  such 
as  those  proposed  in  the  GlueX  project,  are  necessary  to  clarify  this  situation. 


102 


A -9 

o3  3 

-  o 

.2  8 

M  5 

a  ^ 

C3  -4J 


"-s 


g 


;/: 
00 


CO 


u 
$ 

H 
S-< 

o  -o 


a3 


O 

u 

O 
x 


,     CP 

ft   w    C 
a;    o 


o3         t— h 


8  ? 


co     — i   >f-i 


03 

CO 


tf-8 

.5   a 
.3    3  .2 


93 


H  3 


ol 

a)    (B   -p 
.  EH  =£ 

CO 


4- 

— 


bC 


a 

co 

z 


n      03      .. 

"    o    (j 

3  u  a 

o3    a?  •-- 

~g    3    03 

Mi 


DC 
X 

o 


">-<      CD 


a  ^  in 

-  a  P 

a  3  . 

w  53  3 

00  -3  ^ 

«  s  § 

■3  '-S  §: 

H  ■£  < 


ao 

o 

a 

03 

K       ^ 

b 

b 

b 

ca 

CO     ,-H 

e© 

H    CM 

— H 

co 

i     i     i 

'3 

"*    A 

CM 

A  rf 

oi 

^* 

b0 

•  i— 1 

V2 

CD 

lO 

!"H 

O 

410 

o 

o 

O 

otcSS 
2  op  q 

q 

r-sl 

© 

© 

idoo 

O 

* — l*"^ 

i    -H 

-H 

f  1  +  1 

i 

+1 

i     i     > 

Hi 

CT> 

OJ 

°  s 

T-H 

o 

H 

°d 

t-1 

© 

© 

q 

© 

o 



1 

1 

o 

t— 1 

^ 

T—i 

o 

X 

q 

CM 

q 

T— H 

i—i 

_ 

I-H    © 

© 

Oi    o 

O 

o 

.2 

o  -H 

-H 

o  -H 

-H 

-H 

oo  co  i— • 

*-u 

-H  q 

CO 

-H  y3 

^* 

co 

O     ©     r-H 

o 

c3 

»t-  T-i 

DOS  ""I 

T— i 

© 

©     O     © 

M 

°°  co 

^  co 

T— 1 

rH 

h  d  O 

o 

© 

ft 

bO 
3 

o 

-H 

co  -H 

rf    00 

41 

-H 

i—i 

V   V   V 

co 

CO 

o 

i-H 

©' 

s 

03 

M 

n 

1 

+ 

1 
+ 

+   +     a, 

*%  T 

o"^— .  CM 

'a; 
□ 

o3 

o 

T 
38  cm 

i.O 

o 

CM 

< 
<- 

T 

F 

O 
CO 

1— i 

1 — i 

CM    O    CM 
£5    CM    CM 

CM    CM      TT 

"~5   * ^^ 

f- 

.  o 

APPENDIX 
HELICITY  FORMALISM  FOR  TWO  BODY  DECAYS 

Helicity  formalism  can  be  used  to  obtain  the  angular  distribution  of  a  decay 

process.  In  this  appendix  we  restrict  ourselves  to  the  case  of  two  body  decays. 

The  situation  is  as  follows:  particle  A  with  spin  J  a,  and  z-axis  spin  projection  Ma 

decays  to  two  particles,  A  — ►  B  +  C  with  definite  helicity  XB  and  Ac-  Let  the 

direction  of  the  of  the  decay  be  h(8,  (f>)  —  h(9B,  0B)  =  Pb  —  —pc  (see  Figure   1). 

The  amplitude  for  this  decay  is, 

4&a>  < .Wb*o\B\ JaMa  >  ■  (1) 

We  have  little  chance  of  directly  calculating  1.  But  we  can  exploit  the  fact  that 
H  conserves  angular  momentum  by  considering  the  basis  {<  jtuXb^cI}-.  where  j 
is  the  total  angular  momentum  of  B  +  C  and  in  its  projection  along  the  z  axis1  . 
Inserting  this  basis  into  Equation  1, 

AXb\c    =    J2  <  0(t>\BXc\jm\B\c  ><  jm\B\c\H\JAMA  > 

jm 

=    iL,  <  ^^BAc|.?mAsAc  >  SjjASmMAA\B\c 

jm 

=    <  0(})\b\c\JaMa\b\c  >  AXbXc. 


i 


The  rotational  invariance  of  helicities  is  a  necessary  condition  to  construct 
{<     jm\BXc\}  as  a  basis  for  the  two- particle  state  B   +  C.  Rotational  invariance 
is  one  of  two  key  ideas  in  helicity  formalism.  The  other  key  idea,  which  is  boost 
invariance,  is>useful  in  sequential  decays. 
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Figure  1:  A-+B  +  C. 

The  above  equation  contains  two  factorized  terms.  The  first  term  is  indepen- 
dent of  the  interaction  H  and  contains  the  angular  distribution  information  we  are 
interested  in.  The  second  term.  AXa\c,  is  usually  called  the  "helicity  coupling  am- 
plitude1' .  It  is  independent  of  the  angular  distribution,  and  contains  the  physics  of 
the  decay  (in  particular  it  depends  on  \b  and  Ac,  but  not  on  Ma)-  In  the  helicity 
formalism,  the  helicity  coupling  amplitudes  are  unknown  parameters. 

For  strong  and  electromagnetic  decays  there  is  another  invariance  of  H  we 
can  exploit;  parity.  For  these  kind  of  interactions  it  can  be  shown  by  inserting  the 
parity  operator  that, 


A_Ab._Ac  «  (-1) 


Ja—Jb—Jc 


VaVbVcAxbXc 


(2) 


where  7]Ar)sVc  is  the  product  of  the  parity  of  particles  A,  B  and  C.  This  useful 
property  of  A\B\C  reduces  the  number  of  unknown  parameters  in  the  angular 
distribution  of  strong  and  electromagnetic  decays. 

The  next  step  is  to  derive  explicit  formulae  for  <  94>\b\c\  JaM^Xb^c  >  for 
the  B  +  C  system.  Such  derivation,  which  is  out  of  the  scope  of  this  appendix,  will 
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be  omitted  here2  .  The  result  can  be  expressed  in  terms  of  the  Wigner  functions, 

<  6<I>\b\c\JaMa\b\c  >*  ^^^D^^e,-^  (3) 

where  A  =  A#  —  Ac  and, 

The  functions  d^  A(#)  have  real  values,  and  some  of  them  are  given  in  the  particle 
data  book. 

We  have  finally  arrived  to  the  following  important  relation, 

aMa     _    /2  J  a  +  1  D*jA   A  /-x 

All  helicity  distribution  calculations  are  reduced  to  handling  these  basic  ampli- 
tudes. Here  are  some  of  their  very  important  properties, 


K\c  =  °    lf     I^B  -  Ac|  >  7,4,  (6) 


A~ma    ||_.j|  jMA       II  (ys 


The  interested  reader  can  find  this  derivation  in   [35] 
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■  ■ 

Property  6  is  a  reflection  of  the  physical  fact  that  the  projection  of  the  total 
angular  momentum  along  a  particular  axis  cannot  be  greater  than  the  total  angular 
momentum  itself. 

II  ^xB\c  II2  *s  ^e  anSluar  distribution  of  a  final  state  where  MA  and  the 
helicities  of  B  and  C  are  known  (see  Equation  1).  However,  it  is  usually  the  case 
that  MaXbXq  are  not  known.  In  such  cases  we  can  only  measure  the  overall 
angular  distribution,  which  is  simply  the  average  over  possible  initial  states  of  the 
incoherent  sum  of  final  states, 

£-  E  a*^-  E  ^m&JP  w 

MAXBXC  7  MAXBXC 

where  Pma  is  the  weight  of  the  initial  state  Ma  and  dVis  —  d(cos(6B))d(pB- 

Next,  consider  a  sequential  decay,  A  -*  B  4-  C  with  C  -*  D  4-  E  (see  Figure 
2).  The  rotational  and  boost  invariance  of  helicities  makes  it  easy  to  extend  the 
helicity  formalism  to  this  case.  To  do  so  we  consider  a  new  reference  frame,  O', 
at  rest  with  C,  such  that  z'  --  b(#b,$b)3  .  Ail  particle  helicities  remain  the  same 
in  0',  and  in  particular  MZ(C)  --    -\c  (the  negative  sign  arises  from  the  fact 
that  n(#B,  4>b)  is  the  direction  in  which  B  is  traveling).  It  is  now  straightforward 
to  calculate  the  intermediate  amplitude  CX^(E  at  O'.  Of  course  A^B\C  is  still 
calculated  at  ,4's  original  rest  frame.  Since  Ac  is  an  intermediate  state,  the  helicity 
amplitude  for  the  final  state  XbXd^e  is, 


AC 


3  This  condition  does  not  unlikely  determine  0'.  We  can  be  more  precise  by 
defining  O'    -    L(pc)R(<t>B,  &b,  0)0,  where  L  is  the  Lorenz  Boost  operator,  R  is  the 
rotation  operator  in'terms  of  the  Euler  angles,  and  O  is  As  original  rest  frame. 
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Figure  2:  A  —  >  B  ■+  C  and  C  ->  D  +  E.  The  drawing  is  peculiar  in  that  C's  daugh- 
ters are  drawn  as  seen  in  C's  rest  frame  (back -to-back). 

where  we  have  coherently  summed  over  the  intermediate  state  C.  If  Ma^b^d^e 
are  not  known. 


dP 


dtledftn 


c:-    E   a 


MA 


dPx\  a 
dQgdtljj 


£   ^n-^wi 


(12) 


The  sequential  decay  amplitudes  have  the  following  important  properties, 


AC 


-Ma 

AnAnAr 


ACMa 


(13) 


/   #e  /   d*D  ||  a^a,  ||2=  ^  £  ii  <;Ac  ||2||  c%Xb 


Il2 


(14) 


This  kind  of  logic  can  be  extended  to  a  general  two-body  decay  chain. 
Example:  T(IS)  ->  7/2(1270)(/2(1270)  -♦  ir^if) 
From  our  previous  discussion  we  can  immediately  write  (see  Equation    11), 


AC^1S\ 

AT"\r+-\r- 


=   E  A 


Mr(lS)  £-*/2(1270) 

•M/2(i270)    \r+\- 


\f2(1270) 


■•  * 
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In  this  case,  A»-  =  A„.+  =  0,  A7  =  ±1.  In  CLEO  MT(is)  =  ±1*  there  is  no 
Mr(is)  —  0  since  at  the  T(L9)  energy  electrons  couple  only  to  positrons  of  the 
opposite  helicity,  this  way  Pi  —  P_i  =  \  and  Pq  -  0  is  an  excellent  approximation.. 
There  are  therefore  four  amplitudes  we  must  consider,  AC\(j0.  AC{^,  AC[  100  and 


ACT1 


loo- 


The  overall  angular  distribution  is  (see  Equation   12), 


_dP  _ 
di~l1d£Lir+ 


\  II  AC\m  ||?  +1  ||  ACl 


:il  ^ioo  il   +  II  AC_-m 


-100 
II 2 


I!2  +\  I!  AC^  ||2  +~  ||  AC^  f 


where  we  have  used  property   13.  The  two  relevant  amplitudes  are. 


(15) 


AL~m    -—    A^Cqq -\-  AnCm  -f  AUC 


oo 


Arl        —    ax    r®  -a-  a1     rl  a.  a1     p'1 


00 


where  we  have  used  the  fact  that  |A7  —  A/2(i27o)|  <  Jr{is)  ~  1  to  sum  over  all 
possible  A/2{1270). 

All  that  remains  is  to  substitute, 


-4^100  II   : 


A\QC00G  f  -f  j|  A\XC^  ||2  +  |j  -J^fig  f 


f 


AC1       li2 

"HU-100  II 


l-/H10L/00-/tll°00      ^  >H10u00-/H12l-00       r"  -^H^OO  >H12L'00       '    C>C-J 

=  il  ^-ioC  II2  +  II  AU_Ao  ||2  +  II  A\A_2Cl  f  + 

Ml      pO     a*1       p*l    ,      *1      /i0    j*l       /»*2    ,     jl         /?1     4*1       n*1    ,    -_! 
-lO^OO-^-l-l^OO  r  •*  •10l-'00-H-l-  2*-00  ^  -^-1  -  lu00^l-l-2L'00  "+"  '--C-J 


(16) 


into  Equation    15. 

Doing  so  will  give  us  the  full  angular  distribution  of  the  decay  chain.  However, 
using  property   8,  all  the  cross  terms  inside  {...}  vanish  if  we  integrate  over  (f)n+  or 
07.  For  the  time  being  we  integrate  over  both  <j>v+  and  </>7.  We  will  return  to  the 
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terms  inside  {...}  later. 


c2tt  /-2tt  j  p 

Wr+TFTlFr-     -     ^2i(W  Alo  f  +  II  -Alio  f)  II  CSo  IP  + 


/  #>     / 

Jo  Jo 


flKXrd£27r+ 


(II  -ft  ||2  +  ||  A\._x  ||2)  ||  C00  ||2  + 

(il  A II2  +  II  Aw  ii2)  ii  c^  f] 

-    f  [II  AM  ||2|!  Coo  II2  (d\M2  +  4-i{^?)4o{^?  t 
II  An  |j2|i  C00  II2  (AC^)*  +  4>W')4(*r*)1  + 
II  A12  ||2||  Coo  II2  (dU(^)2  +  4-iW2)4>(M2] 

*    f  II  Coo  II2  0!  A10  ||2  \(l  +  cos2  <?7)(|  cos2  0+  -  \f  % 
||  Au  ||2  sin2  07(  J-  sin^+  CO80„+)2  + 
l|A12||2~(l  +  cos2^.)(^sin4^+)]. 

At  this  point  it  is  convenient  to  define  the  normalized  helicity  amplitudes, 

Aio 
a-,Q 


v4!A7o"!l2~H|An||2  +  |[A12|P 

An 

VJtAM^  +  ljAnP  +  ffAjaJP 

Ai2 


G|2 


/||  A10  p  +  I!  Au  P  +  ||  Al2  f 


which  satisfy. 


aio  II2  +  II  an  ||2  +  II  au  ||2=  1.  (17) 


The  normalized  probability  distribution,  ^cos9  w(Cose — )'  can  ^e  0Dtaine^  by 
dividing  the  probability  distribution  by  the  branching  fraction  of  the  sequential 
decay, 
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P    =    2  ||  Coo  II2  (14  IP  +  I!  An  ||2  +  ||Aia  ||2) 

=    I!  Coo  ||2!l  Al0  ||2  +  ||  Coo  ||2||  Au  f  +  \\  Coo  f\\  4s  ||2  + 
II  Coo  !!2||  A_l0  ||2  +  ||  Coo  p||  ii-i-i  ||2  +  ||  Coo  ||2|l  ,4-i-2  ||2 


This  way, 


dU  "  '12  3n    ,    _2fl  Ao 


aio  II    o  (1  +  cos2  ft,) -(3  cos2  ^+  -  1)'  + 


d(cos»7)d(cos0w+)  "        "    8 '  7/8 

3  15 

||  an  ||2  -  sin2  07-—(sin0„.+  cos97r+)2  + 
1  l    4  4 

3  15 

||0i2  ||2  -{I +C0826J-~rW4  fa- 
it is  useful  to  project  the  angular  distribution  on  07  and  0„-+ , 

"  aw  ||2  -(3cos207r+-l)2+  ||  an  ||2  —  (sin0T+  cos07r+)2+  ||  ai2  ||2  — sin4^ 


dcos^f      "    iU"    8V  *       '       "    ""     4V       "  "'       "    ""    16 

(18) 


"^  '"  i|2    ,    II  „       Il2\^/i     , 2/i   \  ,    II  „       ||2  ^    .    o 


=  (II  a10  i|2  +  ||  a12  |i2)^(l  +  cos207)+  ||  an  ||2  Jsin207.  (19) 

a  cos  f-y  o  4 

Let  us  now  return  to  the  terms  inside  {...}  in  Equation! 6.  It  turns  out  that  for 
these  terms  a  substantial  simplification  occurs  if  we  integrate  over  cos  0w->- .  After 
some  algebra,  we  arrive  at  the  following  expression, 


'1  -  ^={ai0a*2  +  a*lQal2}  cos[2(</>7  -  (pn+)}).  (20) 


d(07  -  </v+)      2tt         a/6 

Example:  T(l,9)  -►  7/4(2050)(/4(2050)  ->  tt+tt") 

For  the  case  /4(2050)  the  calculation  is  very  similar  to  the  /2(1270)  case 
discussed  above.  The  allowed  helicities  for  ^(2050)  are  the  same  as  the  allowed 
helicities  for  /2(1270).  This  makes  the  angular  distribution  calculation  very  similar 


Ill 


to  the  /2(1270)  case,  the  only  difference  being  that  the  functions  d£    (1270)o  Me 
replaced  bv  d\  n. 

The  photon's  angular  distribution  comes  out  to  be  the  same.  The  result  for 
the  pion's  normalized  angular  distribution  is. 

"  aio  il2  ~(35cos4^+  -  30  cos2  ^+  +  3)2  + 


dcos6n+  128 

AC 

i|  an  ||2  ^(sin.(?,l+  cos^+)2(7cos2^+  -  3)2  + 

(21) 

Fit  of  the  Helicity  Angular  Distribution  Obtained  from  Data 

in  practice,  a  particular  bin  in  the  helicity  angular  distribution  obtained  from 
data  has  bins  with  dN  =  NdU  events,  where  N  is  the  total  number  of  events 
(which  is  allowed  to  float  during  the  fit),  and  dU  is  the  portion  of  the  normalized 
probability  assigned  to  the  bin  by  the  helicity  formalism  prediction. 

To  incorporate  Equation  17  directly  into  the  fit  we  make  the  following  change 
of  variables, 

j|  aw  ||=  cos  9  cos  $;  ||  au  [|=  cos 0 sin*;  ||  an  ||=  sinO.  (22) 

Because  of  this  definition,  both  6  and  $  should  be  required  to  be  within  the 
interval  (0,  f )  during  the  fit.  The  new  variables  have  the  convenient,  property  that 
for  J=0,  $-0=0,  and  for  J  =  1,  $=^0.  in  general  there  are  1,  2  and  3  degrees 
of  freedom  for  the  fits  with  J  =  0,  J  -  1,  and  J  >  1  respectively.  One  degree 
of  freedom  always  corresponds  to  the  number  of  events,  N,  and  the  other  two 
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correspond  to  0  and  $,  which  determine  the  magnitude  of  the  normalized  helicity 
amplitudes  through  Equation  22  4  . 

When  fitting,  it  is  preferable  to  choose  the  full  2D  helicity  angular  distribution 
which  is  a  function  of  #7  and  8V+  instead  of  its  one  dimensional  projections.  The 
2D  choice  is  better  not  only  because  the  2D  distribution  carries  the  maximum 
amount  of  information,  but  because  in  general  the  efficiency  depends  on  both  the 
pion  and  photon  helicity  angles,  which  can  lead  to  complications  when  efficiency 
correcting  the  one-dimensional  projections.  Building  on  the  two  previous  examples 
the  2D  helicity  angular  distributions  for  ,/  =  0,  1,  2,  3,  4  are  (for  compactness  we 
define  x  —  cos  8V+  and  y  ~  cos#7), 


dNJ=0  3  o,l 


^^W^'j  (23i 


sm'e^l-^U-z2)) 


(24) 


4  Obtaining  the  error  of  the  normalized  helicity  amplitudes  requires  some  matrix 

(ai0\  _       /0\ 

multiplication.  For  compactness  let  us  define  a  —      012     and  cD  = ■  (  ,  1 .  The  error 

\an/  ^     ' 

matrix  for  the  normalized  helicity  amplitudes  is, 

<  S'dSa1  >=  t—  <  5u6u)T  >  7—; 
O'jJ  ou 

where  <  Sio6uJr  >  is  the  error  matrix  obtained  from  the  fit. 


_ 
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^  =  ^2(cos20cos2$  |(1  f  y2)^(3z2  -  1)2+ 

sin26  |(1  -  y2)^!  -  a2)*  (25) 

cos^sin^la  +  ^Cl-a;2)2) 

8  '  lo 


dpfl  =  A-J=s(  cos2  0 cos2  $  ? (1  +  W2)^5a:3  ~  3s)2+ 
dxdy  8  8 

sin2  6  hi-  y2)^(l  -  x2)(5x2  -  1)2+  (26) 

cos2esin2$?(l  +  y2j~x2(l-x2)2) 
8  16 


^%-4  =  JVJ=4(cos2  6cos2  $  |(1  +  J/2)  r^(35x4  -  30x2  +  3)2+ 
dxdy  8  128 

sin2  6  |(1  -  y2)%x2{\  -  x2)(7x2  -  3)2+  (27) 

cos20sin2$  |(l  +  |^)g(l-^)»(7«»  -l)2). 

In  this  analysis,  because  of  our  low  statistics  (and  by  low  statistics  we  mean 
that  there  are  numerous  bins  with  <  1  events  in  the  2D  plots),  we  project  Equa- 
tions 23-27  on  x  =  cos  6^,  and  y  =  cos07.  We  also  fold  each  distribution  axound  its 
symmetry  axis  x  =  0  and  y  —-  0  to  increase  each  bin's  statistics.  Finally,  we  do  a 
simultaneous  fit  to  both  distributions  and  obtain  N,  0,  and  $. 

Admittedly,  this  procedure  has  the  problems  described  previously  when  2D  fits 
where  advocated,  which  can  give  rise  to  systematic  effects  in  the  measured  helicity 
amplitudes.  In  essence,  we  need  to  use  a  unitary  weight  function  proportional 
to  the  efficiency  distribution  (see  Figure   7-3)  in  the  integrals  used  to  make  the 
projections.  This  affects  the  relative  strength  of  each  helicity  amplitude  and  the 
overall  efficiency,  but  does  not  change  each  projected  shape.  However,  the  helicity 
amplitudes  themselves  are  only  used  to  calculate  systematic  effects  on  the  efficiency 
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of  resonances,  so  by  simultaneously  fitting  the  unweighted  ID  projections  we  are 
just  ignoring  the  systematics  on  the  systematica. 

Finally,  let  us  consider  the  07  -  0W+  distribution  for  J  =  2.  Equation  20 
becomes, 


dNj=2  Nj^_2 , 


i 


~ {cos2  6  sin  2$  cos  A}  cos[2(07  -  0,.+  )])        (28) 


rf(<£7  —  (f>n+ )  ITk  \J§ 

where  the  new  degree  of  freedom  A  is  the  relative  phase  between  aio  and  a2o-  The 
deviation  from  the  flat  distribution  is  characterized  by  R  =  -^{cos2  0  sin  2$  cos  A}  < 
4? .  The  bigger  this  factor  is,  the  easier  it  is  to  observe  in  data.  Figure  3  shows  the 
continuum  subtracted  <£7  -  0„+  distribution  for  events  near  the  /2(1270)  resonance. 
The  factor  R  obtained  from  fitting  the  distribution  is  non-significant.  It  is  a  bit 
tedious  but  straightforward  to  calculate  the  <£7  -  (pn+  distribution  for  other  J  values. 
We  will  not  present  those  results  here.  They  are  not  very  interesting  because  we 
need  large  statistics  to  observe  deviations  from  the  flat  distribution,  and  there  is  no 
effect  on  the  efficiency  from  non-fiat  c/>7  —  </>„.+  distributions. 
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MINUIT  x'  Fit  to  Plot  9&3 

U1S->Gamma  h  h.  pipiphi  axis 

File:  Generated  internally 

Plot  Area  Total/Fit  779.39  /  779.39 

Func  Area  Total/Fit  769.52  /  769.52 

y2=     5.2  for  10-  2d.o.f., 

Errors  Parabolic 

Function  1:COMIS  Function  XMNCMJ 
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R  -6.87430E-02  +   0.106? 
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Figure  3:  Continuum  subtracted  ri>7  —  ^+  distribution  for  events  near  the  /2(1270) 
resonance  with  a  fit  to  Equation  28  overkyed. 
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Example:  Continuum  Events. 

Our  previous  examples  tackled  radiative  decays  of  the  T(15).  Now  let  us 
consider  the  continuum  "radiative  beam"  process  e+e~  -*  77r+7r~  (see  Figure  4). 
At  first  glance  it  might  seem  that  the  helicity  formalism  cannot  be  applied  in  this 
case.  But  a  closer  look  at  the  diagram  of  this  process,  shown  in  Figure  4,  indicates 
otherwise.  The  key  point  is  that,  regardless  of  the  particular  coupling,  helicity 
formalism  can  still  be  applied  at  the  CM  of  the  7r+7T'  vertex,  where  the  initial  state 
has  Jp  =  l~. 


7 

Figure  4:  e+e~       — *•       jtt+tx"  in  continuum.  7r4'7r_  pair  production  usually  occurs 
through  p  —  7  coupling. 


We  only  need  to  know  that  we  are  dealing  with  1     — *  0  0    at  the  7r+7r 
vertex.  There  is  only  one  helicity  amplitude  we  need  to  consider, 

■Ajo  =  y^54»  (29) 

which  after  a  couple  simple  steps  leads  to, 


=  iV^sin2^+.  .    (30) 


d  cos  #„.+  4 

This  result  is  also  valid  for  the  kaon  case. 

On  the  other  hand,  the  photon's  angular  distribution  cannot  be  calculated 
using  two-body  helicity  formalism.  A  good  approximation  is  (see  [36]  Appendix  B), 
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dP  l  +  cos207  ,01. 


d  cos  #7       1  —  cos2  #7 ' 

which  is  not  valid  at  low  8  angles  (6  — >  0,  61  — ►  n)  because  of  the  effects  on  the 
approximation  of  the  neglected  electron  mass. 

Figure   5  shows  these  angular  distributions  in  4S  "radiative  beam"  events.  The 
agreement  between  theory  and  experiment  is  excellent.  This  is  reassuring,  because 
in  the  main  analysis  we  rely  on  our  angular  distribution  measurements  to  get  the 
correct  efficiency  in  data  and  to  identify  the  J  value  of  each  resonance. 


\ 


118 


> 

1500 

u 

s 

1000 

1/1 

— 

-^ 

500 

BO 

1 

o 

u- 

7 

c 

0 

0.100 

- 

£ 

'fi 

'X. 

0.050 

0 

> 

0.000 

lil 

= 

0 

0.30 

3 

Jl 

'~ 

0.20 

t» 

a 

0.10 

= 

> 

0.00 

u 


0.0 


7t  7i"  Invariant  mass  (GeV/c  ) 


K+K  Invariant  mass  (GeV/c  ) 


Figure  5:  Helicity  angle  distribution  for  "radiative  beam"  continuum  events 
e+e"     -->    7p  and  e+e~     —>    j(f>  in  the  4S  data.  Top  left  (right)  plot  shows  in  gray 
the  events  selected  as  p    —*    7r+7r~  (4>    — »    K+K~).  Middle  left  (right)  plot  shows 
the  7r+7r~  (K+K~)  helicity  angle  distribution  for  the  selected  events,  and  has  the 
helicity  formalism  prediction,  Equation  30,  overlaid.  Bottom  left  (right)  plot  shows 
the  7  angular  distribution  for  the  photon  accompanying  the  p  (</>),  and  has  the  first 
order  QED  prediction,  Equation  31,  overlayed. 
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